00 

9 

cr 


DC 

LU 

O 


o 

Q 

DC 

LU 


US  Army  Corps 
of  Engineers® 

Engineer  Research  and 
Development  Center 


Source  Term  Model  for  Fine  Particle 
Resuspension  from  Indoor  Surfaces 


Yoojeong  Kim,  Ashok  Gidwani,  Mark  Sippola,  and  Chang  W.  Sohn  February  2008 


tuo 

c 

■ 

0  >, 

c  o 

c*  2 

LU  O 
O  _| 


o 


o 


+-»  CO 
(/)  0 
C  0 
O  0 

o  cr 


Approved  for  public  release;  distribution  is  unlimited. 


ERDC/CERL  TR-08-4 
February  2008 


Source  Term  Model  for  Fine  Particle 
Resuspension  from  Indoor  Surfaces 


Yoojeong  Kim 
Triton  Systems,  Inc. 

200  Turnpike  Road 
Chelmsford,  MA  01824 

Ashok  Gidwani 
CFD  Research  Corporation 
215  Wynn  Drive 
Huntsville,  AL  35805 

Mark  Sippola 

Lawrence  Berkeley  National  Laboratory 
1  Cyclotron  Road 
Berkeley,  CA  94720 

Chang W.  Sohn 

U.S.  Army  Engineer  Research  and  Development  Center 
Construction  Engineering  Research  Laboratory  (CERL) 
2902  Newmark  Dr. 

Champaign,  IL  61824 

Final  Report 


Approved  for  public  release;  distribution  is  unlimited. 


Prepared  for  U.S.  Army  Corps  of  Engineers 
Washington,  DC  20314-1000 

Under  Work  Units  F726C1  and  J5LG71 


ERDC/CERL  TR-08-4 


n 


Abstract:  This  Phase  I  effort  developed  a  source  term  model  for  particle 
resuspension  from  indoor  surfaces  to  be  used  as  a  source  term  boundary 
condition  for  CFD  simulation  of  particle  transport  and  dispersion  in  a 
building.  Specifically,  this  work:  (l)  investigated  responsible  mechanisms 
for  fine  particle  resuspension  from  indoor  surfaces,  (2)  identified 
parameters  relevant  to  resuspension,  (3)  performed  a  dimensional 
analysis  and  derivation  of  a  resuspension  model,  and  (4)  evaluated  the 
model  against  published  experimental  data  on  resuspension.  Preliminary 
validation  of  the  derived  model  was  conducted  based  on  a  set  of 
experimental  data  from  the  Lawrence  Berkeley  National  Laboratory. 


DISCLAIMER:  The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or  promotional  purposes.  Citation 
of  trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of  such  commercial  products.  All  product 
names  and  trademarks  cited  are  the  property  of  their  respective  owners.  The  findings  of  this  report  are  not  to  be  construed  as 
an  official  Department  of  the  Army  position  unless  so  designated  by  other  authorized  documents. 

DESTROY  THIS  REPORT  WHEN  NO  LONGER  NEEDED.  DO  NOT  RETURN  IT  TO  THE  ORIGINATOR. 


ERDC/CERL  TR-08-4 


Contents 

Figures  and  Tables . v 

Preface . viii 

Unit  Conversion  Factors . ix 

1  Introduction . 1 

1.1  Background . 1 

1.2  Objective . 2 

1.3  Approaches . 2 

1.4  Scope . 2 

1.5  Mode  of  Technology  Transfer . 3 

2  Resuspension  of  Particles . 4 

2.1  Resuspension  Mechanism . 5 

2.1.1  Adhesion  Forces  5 

2.1.2  Removal  Forces  7 

2.2  Existing  Models  of  Resuspension . 9 

2.2.1  Resuspension  Factor  K  and  Resuspension  Rate  A  9 

2.2.2  Empirical  Models  9 

2.2.3  Theoretical  Models  11 

2.3  Responsible  Parameters  for  Resuspension . 13 

3  Development  of  Indoor  Resuspension  Models  of  Fine  Particles . 15 

3.1  Identification  of  Relevant  Parameters . 15 

3.2  Dimensional  Analysis . 16 

3.2.1  Case  1  16 

3.2.2  Case  2  16 

3.2.3  Case  3  17 

3.2.4  Case  4  17 

3.3  Model  Development . 17 

3.3.1  Empirical  Correlations  18 

3.3.2  Physics-Based  Model  27 

3.3.3  Comparison  of  Dimensionless  Variables  in  Empirical  Correlations  and  Physics- 

based  Model  30 

3.4  Testing  of  Models . 31 

3.4.1  Collected  Experimental  Data  31 

3.4.2  Model  Testing  34 

4  Source  Term  Model  for  CFD  Simulation . 38 

4.1  Description  of  Particle  Resuspension  in  a  Bench-top  Test  Chamber . 38 

4.2  CFD  Simulations  of  Flow  Field  in  the  Test  Chamber . 39 

4.2.1  Experimental  Results  39 

4.2.2  CFD  Simulations  of  Flow  Field  40 


ERDC/CERL  TR-08-4 


IV 


4.3  CFD  Simulation  of  Particle  Resuspension  and  Transport . 47 

4.3.1  Resuspension  Rate  47 

4.3.2  CFD  Simulations  of  Particle  Transport  and  Resuspension  48 

4.4  Comparisons  of  Resuspension  Data  between  Simulations  and  Measurement . 49 

5  Summary,  Conclusions,  and  Recommendations . 58 

5.1  Summary . 58 

5.2  Conclusions . 58 

5.3  Recommendations . 59 

References . 60 

Acronyms,  Abbreviations,  and  Terms . 63 

Report  Documentation  Page . 66 


ERDC/CERL  TR-08-4 


v 


Figures  and  Tables 

Figures 

1  Gantt  chart  detailing  Phase  I  study  tasks . 2 

2  Van  der  Waals  force  between  particle  and  plane . 5 

3  Van  der  Waals  force  between  planes . 6 

4  Resuspension  rate  versus  time  plotted  from  the  data  of  Ibrahim  (2004)  and 
Nicholson  (1993).  For  the  corresponding  time,  (a)  used  the  mid-point  and  (b) 

used  t  defined  by  Eq  45 . 20 

5  Plot  of  \n(AdP/u*)  vs.  In (u*t/dP)  of  all  SS70  results  with  different  RH  from 

Ibrahim  (2004) . 23 

6  Resuspension  rate  versus  time  plotted  from  the  data  used  to  test  empirical 
correlations.  For  the  corresponding  time,  (A)  used  the  mid-point  and  (B)  used 

f  defined  by  Eq  45 . 32 

7  Performance  of  the  empirical  correlations,  (A)  EC  I,  (B)  EC  IV,  and  (C)  EC  VII, 

against  experimental  data . 35 

8  A  schematic  diagram  of  the  chamber  for  particle  resuspension  experiments; 
all  dimensions  are  in  inches,  and  the  x-y  positions  of  the  measurement  points 

for  velocity . 39 

9  (A)  3D  model  of  chamber  showing  the  general  direction  of  air  flow  through 
various  inlets  and  outlets.  (B)  unstructured  (binary-tree)  mesh  implemented 
inside  the  chamber  3D  model.  The  refined  mesh  layer  just  below  the  make-up 

air  vent  is  to  resolve  the  flow  modification  through  the  two  fans . 41 

10  Snapshot  of  the  CFD-ACE+  graphical  user  interface  with  physical  models 

used  in  the  simulations . 42 

11  (A)  Velocity  vectors  colored  by  magnitude  of  velocity  in  the  plane  through  the 
fan  centers;  vectors  are  also  plotted  in  the  plane  along  the  diagonal  through 
bottom;  isosurfaces  of  velocity  magnitude  in  the  plane  of  the  two  fans 

(6  L/min);  (B)  isocontours  of  the  friction  velocity  at  the  floor  (6  L/min) . 44 

12  (A)  Velocity  vectors  colored  by  magnitude  of  velocity  in  the  plane  through  the 
fan  centers;  vectors  are  also  plotted  in  the  plane  along  the  diagonal  through 
bottom;  isosurfaces  of  velocity  magnitude  in  the  plane  of  the  two  fans 

(20  L/min);  (B)  isocontours  of  the  friction  velocity  at  the  floor  (20  L/min) . 45 

13  (A)  Velocity  vectors  colored  by  magnitude  of  velocity  in  the  plane  through  the 
fan  centers;  vectors  are  also  plotted  in  the  plane  along  the  diagonal  through 
bottom;  isosurfaces  of  velocity  magnitude  in  the  plane  of  the  two  fans 

(40  L/min);  (B)  isocontours  of  the  friction  velocity  at  the  floor  (40  L/min) . 45 

14  (A)  Velocity  vectors  colored  by  magnitude  of  velocity  in  the  plane  through  the 
fan  centers;  vectors  are  also  plotted  in  the  plane  along  the  diagonal  through 
bottom;  isosurfaces  of  velocity  magnitude  in  the  plane  of  the  two  fans 

(60  L/min);  (B)  isocontours  of  the  friction  velocity  at  the  floor  (60  L/min) . 46 

15  (A)  Velocity  vectors  colored  by  magnitude  of  velocity  in  the  plane  through  the 
fan  centers;  vectors  are  also  plotted  in  the  plane  along  the  diagonal  through 


ERDC/CERL  TR-08-4 


VI 


bottom;  isosurfaces  of  velocity  magnitude  in  the  plane  of  the  two  fans 

(132  L/min);  (B)  isocontours  of  the  friction  velocity  at  the  floor  (132  L/min) . 46 

16  Particle  source  from  the  floor  as  a  function  of  time.  The  particle  source  was 
calculated  using  the  resuspension  rate  given  by  the  empirical  correlation,  EC 
I.  The  particle  size  was  0.5  jum,  and  the  initial  deposition  amount  was  352  ng 

(0.06  ng/cm2).  The  floor  material  was  linoleum . 49 

17  Isocontours  of  particle  concentration  in  a  plane  at  the  floor  and  in  a  plane 
along  the  diagonal  of  the  chamber  showing  spatial  and  temporal  evolution  of 
particle  concentration  between  4  and  30  seconds:  at  20  L/min  (left)  and 
132  L/min  (right).  The  particle  size  was  0.5  jum,  and  the  initial  deposition 

amount  was  352  ng  (0.06  ng/cm2).  The  floor  material  was  linoleum . 50 

18  Isocontours  of  particle  concentration  in  a  plane  at  the  floor  and  in  a  plane 
along  the  diagonal  of  the  chamber  showing  spatial  and  temporal  evolution  of 
particle  concentration  between  60  and  240  seconds:  at  20  L/min  (left)  and 
132  L/min  (right);  the  particle  size  was  0.5  jum,  and  the  initial  deposition 

amount  was  352  ng  (0.06  ng/cm2);  the  floor  material  was  linoleum . 51 

19  Time  evolution  of  the  particle  concentration  at  the  outlet  located  5.25  in. 
above  the  floor  at  the  center  between  (A)  0  -  600  s  and  (B)  0  -  3600  s;  the 
particle  size  was  0.5  |um,  and  the  initial  deposition  amount  was  352  ng  (0.06 


ng/cm2);  the  floor  material  was  linoleum . 52 

20  The  initial  mass  deposition  density  as  a  function  of  particle  size  on  (A) 

linoleum  and  (B)  carpet . 54 


21  The  total  fraction  of  mass  resuspended  during  1  hr  of  resuspension 
experiment  for  (A)  linoleum  and  (B)  carpet.  The  fraction  of  mass  resuspended 
was  calculated  for  a  given  particle  size  as  the  ratio  of  the  mass  of  particle 
resuspension  to  the  mass  of  initial  deposition  on  the  floor.  The  fan  setting 
was  75,  and  the  measured  fan  flow  rate  was  20  L/min.  The  red  and  the  blue 
circles  indicate  predictions  using  EC  I  at  20  L/min  and  132  L/min, 
respectively,  of  the  fan  flow  rate.  The  purple  circles  represent  the  prediction 
using  the  empirical  correlation  given  by  Loosmore  (2003)  using  20  L/min  as 

the  fan  flow  rate . 54 

22  The  total  fraction  of  mass  resuspended  during  1  hr  of  resuspension 
experiment  for  (A)  linoleum  and  (B)  carpet.  The  fraction  of  mass  resuspended 
was  calculated  for  a  given  particle  size  as  the  ratio  of  the  mass  of  particle 
resuspension  to  the  mass  of  initial  deposition  on  the  floor.  The  fan  setting 
was  110,  and  the  measured  fan  flow  rate  was  40  L/min.  The  red  and  the  blue 
circles  indicate  predictions  using  EC  I  at  40  L/min  fan  flow  rate,  and  tstart  and 
tend  specify  the  lower  limit  and  the  upper  limit,  respectively,  of  the  integration 

in  Eq  80 . 54 

23  The  total  fraction  of  mass  resuspended  during  1  hr  of  resuspension 
experiment  for  (A)  linoleum  and  (B)  carpet.  The  fraction  of  mass  resuspended 
was  calculated  for  a  given  particle  size  as  the  ratio  of  the  mass  of  particle 
resuspension  to  the  mass  of  initial  deposition  on  the  floor.  The  fan  setting 
was  140,  and  the  measured  fan  flow  rate  was  70  L/min.  The  red,  the  blue, 
and  the  purple  circles  indicate  predictions  using  EC  I  at  60  L/min  of  the  fan 
flow  rate,  and  tstart  and  tend  specify  the  lower  limit  and  the  upper  limit, 
respectively,  of  the  integration  in  Eq  80 


55 


ERDC/CERL  TR-08-4  vii 


Tables 

1  Models  for  the  resuspension  factor  K . 10 

2  Experimental  parameters  for  works  by  Ibrahim  (2004)  and  Nicholson  (1993) . 19 

3  Hamaker  constants  used  for  calculating  A132  for  Ibrahim’s  and  Nicholson’s 

data . 21 

4  Variables  for  the  experimental  work  by  Ibrahim  (2004) . 22 

5  The  multiple  regression  result  of  Eq  51  using  Ibrahim’s  (2004)  and 

Nicholson’s  (1993)  data  with  the  mid-point  as  the  corresponding  time  for  the 
resuspension  rate . 24 

6  Empirical  correlations  obtained  using  different  groups  of  dimensionless 
numbers.  Each  dimensionless  function  was  evaluated  using  both  the  mid¬ 
point  and  f  for  the  time  that  represents  the  resuspension  rate,  and  the  effect 

of  the  RH  was  assumed  to  be  negligible . 26 

7  Dependency  of  A  on  each  variable  for  different  models . 26 

8  Experimental  parameters  for  works  by  Giess  (1997)  and  Braaten  (1990) . 32 

9  Experimental  parameters  for  works  by  Wu  (1992),  Krauter  (2007),  and  Wen 

(1989) . 33 

10  Additional  Hamaker  constants  used  for  calculating  A132  for  Giess  et  al. 

(1997),  Braaten  et  al.  (1990),  Wu  et  al.  (1992),  Krauter  and  Biermann  (2007), 

and  Wen  and  Kasper  (1989) . 34 

11  Sensitivity  analyses  of  the  empirical  correlations  of  EC  I,  EC  IV,  and  EC  VII . 37 

12  Measured  flow  rates  of  the  variable  speed  fan . 40 

13  The  average  of  the  minimum,  the  average,  and  the  maximum  air  velocity  in 
the  horizontal  direction  measured  at  the  nine  measurement  points  in  the 

plane  10  cm  above  the  floor  for  the  four  different  fan  speeds . 40 

14  The  average  of  the  minimum,  the  average,  and  the  maximum  air  velocity  in 
the  vertical  direction  measured  at  the  nine  measurement  points  in  the  plane 

10  cm  above  the  floor  for  the  four  different  fan  speeds . 40 

15  Measured  and  calculated  fan  flow  rates;  the  average  fan  flow  velocities  were 
calculated  from  the  measured  fan  flow  rates  from  the  experiments  using  the 

outer  radius  and  the  hub  radius  of  the  fans . 43 

16  Measured  and  calculated  average  velocities  of  the  nine  measurement  points 

in  the  plane  10  cm  above  the  floor  for  various  fan  flow  rates;  the  average 
horizontal  and  vertical  velocities  are  from  Tables  13  and  14 . 43 

17  Measured  and  calculated  average  velocities  of  the  nine  points  in  the  plane 
that  the  fans  were  placed  (10  in.  above  the  floor)  for  various  fan  flow  rates; 
the  x-y  positions  of  the  nine  points  were  the  same  as  those  of  the  nine 
measurement  points  in  Tbl.  16,  but  the  z-positions  for  these  points  were  10-in 

from  the  floor  instead  of  10  cm . 43 

18  The  homogenous  Hamaker  constants  of  the  floor  materials  and  the  particle 

used  for  calculating  A132  for  the  chamber  experiments . 47 

19  Properties  of  particles  as  a  function  of  the  particle  size 


48 


ERDC/CERL  TR-08-4 


viii 


Preface 

This  study  was  conducted  for  the  Engineer  Research  and  Development 
Center,  Construction  Engineering  Research  Laboratory  (ERDC-CERL) 
under  MIPRNo.  50633728  (FY05),  “Experimental  Determination  of  Par¬ 
ticle  Resuspension  from  Room  Surfaces  Due  to  Air  Turbulence”  and  under 
Phase  I  SBIR  Contract  No.  W9132T-07-C-0013,  “Source  Term  Model  for 
Fine  Particle  Resuspension  from  Indoor  Surfaces,”  Work  Unit  F726C1, 
“Source  Term  Model  for  Fine  Particles.”  The  technical  Monitor  was  Martin 
J.  Savoie,  CEERD-TD. 

The  work  was  performed  by  the  Energy  Branch  (CF-E)  of  the  Facilities  Di¬ 
vision  (CF),  Construction  Engineering  Research  Laboratory  (CERL).  The 
CERL  Principal  Investigator  was  Dr.  Chang  W.  Sohn.  This  work  was  partly 
done  by  the  Lawrence  Berkeley  National  Laboratory,  Berkeley,  CA  and  Tri¬ 
ton  Systems,  Inc.,  Chelmsford,  MA.  Dr.  Thomas  A.  Hartranft  is  Chief, 
CEERD-CF-E,  and  L.  Michael  Golish  is  Chief,  CEERD-CF.  The  associated 
Technical  Director  was  Mr.  Martin  J.  Savoie,  CEERD-TD.  The  Director  of 
CERL  is  Dr.  Ilker  R.  Adiguzel. 

CERL  is  an  element  of  the  U.S.  Army  Engineer  Research  and  Development 
Center  (ERDC),  U.S.  Army  Corps  of  Engineers.  The  Commander  and  Ex¬ 
ecutive  Director  of  ERDC  is  COL  Richard  B.  Jenkins,  and  the  Director  of 
ERDC  is  Dr.  James  R.  Houston. 


ERDC/CERL  TR-08-4 


IX 


Unit  Conversion  Factors 


Multiply 

By 

To  Obtain 

cubic  feet 

0.02831685 

cubic  meters 

cubic  inches 

1.6387064  E-05 

cubic  meters 

cubic  yards 

0.7645549 

cubic  meters 

degrees  (angle) 

0.01745329 

radians 

degrees  Fahrenheit 

(F-32)/1.8 

degrees  Celsius 

feet 

0.3048 

meters 

gallons  (U.S.  liquid) 

3.785412  E-03 

cubic  meters 

inches 

0.0254 

meters 

inch-pounds  (force) 

0.1129848 

newton  meters 

miles  (nautical) 

1,852 

meters 

miles  (U.S.  statute) 

1,609.347 

meters 

miles  per  hour 

0.44704 

meters  per  second 

mils 

0.0254 

millimeters 

ounces  (mass) 

0.02834952 

kilograms 

ounces  (U.S.  fluid) 

2.957353  E-05 

cubic  meters 

pints  (U.S.  liquid) 

4.73176  E-04 

cubic  meters 

pints  (U.S.  liquid) 

0.473176 

liters 

pounds  (force) 

4.448222 

newtons 

pounds  (force)  per  foot 

14.59390 

newtons  per  meter 

pounds  (force)  per  inch 

175.1268 

newtons  per  meter 

pounds  (force)  per  square  foot 

47.88026 

pascals 

pounds  (force)  per  square  inch 

6.894757 

kilopascals 

pounds  (mass) 

0.45359237 

kilograms 

quarts  (U.S.  liquid) 

9.463529  E-04 

cubic  meters 

square  feet 

0.09290304 

square  meters 

square  inches 

6.4516  E-04 

square  meters 

square  miles 

2.589998  E+06 

square  meters 

square  yards 

0.8361274 

square  meters 

tons  (2,000  pounds,  mass) 

907.1847 

kilograms 

yards 

0.9144 

meters 

ERDC/CERL  TR-08-4 


1 


1  Introduction 

1.1  Background 

How  contaminants  are  dispersed  inside  buildings  has  long  been  a  topic  of 
interest  to  those  seeking  to  improve  indoor  air  quality  (IAQ)  in  modern 
buildings.  This  topic  also  interests  those  seeking  to  protect  against  the  po¬ 
tential  for  indoor  dispersal  of  toxic  contaminants,  i.e.,  chemical  and  bio¬ 
logical  (CB)  agents. 

To  design  an  active  protection  system  or  to  develop  successful  counter¬ 
measure  strategies  against  CB  agents  dispersed  indoors  requires  detailed 
information  on  a  room’s  dispersion  profile.  A  number  of  computational 
fluid  dynamics  (CFD)  codes  capable  of  simulating  contaminant  dispersion 
inside  a  building  already  exist  in  the  public  and  commercial  domains 
(Gadgil  et  al.  2000;  Sohn  et  al.  2004).  However,  to  apply  these  CFD  codes 
to  calculate  a  dispersion  profile  first  requires  a  prescription  of  boundary 
conditions;  one  of  the  critical  boundary  conditions  affecting  the  simulation 
of  contaminant  dispersion  profile  is  the  “source  term  modeling.” 

In  fiscal  year  2007  (FY07),  the  Army  Small  Business  Innovative  Research 
(SBIR)  program  office  funded  a  study  of  source  term  modeling  for  fine 
particle  resuspension  from  indoor  surfaces.*  This  Phase  I  study: 

1.  Investigated  responsible  mechanisms  for  fine  particle  resuspension 
from  indoor  surfaces 

2.  Identified  parameters  relevant  to  resuspension 

3.  Performed  dimensional  analysis  and  derivation  of  a  resuspension 
model 

4.  Evaluated  the  model  against  published  experimental  data  on  resuspen¬ 
sion. 

This  work,  based  on  the  results  of  the  SBIR  Phase  I  contract  research  per¬ 
formed  from  April  through  October  2007,  was  undertaken  to  develop  a 
source  term  model  for  particle  resuspension  under  indoor  environment. 


*  Chemical  and  Biological  Defense  SBIR  Program  2007.1  Solicitation,  Topic  Number:  CBD07-109, 
Source  Term  Model  for  Fine  Particles  off  Indoor  Surfaces,  accessible  through  URL: 
https://www.armvsbir.com/portal/cbd071/menu.asp 
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1.2  Objective 

The  overall  objective  of  this  Phase  I  effort  was  to  develop  a  source  term 
model  for  particle  resuspension  under  indoor  environment  to  be  used  as  a 
source  term  boundary  condition  for  CFD  simulation  of  particle  transport 
and  dispersion  in  a  building.  The  specific  objectives  were  to: 

1.  Determine  the  responsible  mechanisms  for  fine  particle  resuspension 
from  indoor  surfaces 

2.  Identify  governing  parameters  relevant  to  resuspension  off  indoor  sur¬ 
faces 

3.  Develop  models  for  resuspension  of  fine  particles 

4.  Evaluate  the  models  against  experimental  data  in  the  public  domain. 

1.3  Approaches 

This  work  began  with  an  exhaustive  literature  study  of  published  work  on 
the  related  subject.  Figure  1  shows  the  six  project  tasks  and  the  time  in 
which  they  were  accomplished. 

1.4  Scope 

This  report  describes  work  conducted  for  the  Phase  I  activities  described 
in  Figure  1.  The  validation  of  the  model  against  the  controlled  experiments 
in  a  prototype  test  chamber  is  beyond  the  scope  of  the  Phase  I  work.  The 
validation  of  the  derived  model  through  a  series  of  experiments  in  a  proto¬ 
type  test  chamber  is  included  in  the  Phase  II  SBIR  proposal. 


Tasks 

Months 

1 

2 

3 

4 

5 

6 

Task  1:  Investigate  the  responsible 

mechanisms  for  particle  resuspension 

Task  2:  Identify  relevant  parameters  for 
resuspension 

— 

Task  3:  Perform  dimensional  analysis  and 
derive  an  analytic  model 

■ 

Task  4:  Compile  experimental  data  on  particle 
resuspension 

Task  5:  Compare  the  model  to  experimental 
data 

■ 

Task  6:  Perform  demonstration  CFD  simulation 
with  the  model  as  a  B.C. 

Figure  1.  Gantt  chart  detai 


ing  Phase  I  study  tasks. 
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1.5  Mode  of  Technology  Transfer 

It  is  anticipated  that  the  results  of  this  work  will  be  submitted  to  a  peer- 
reviewed  journal.  This  report  will  be  made  accessible  through  the  World 
Wide  Web  (WWW)  at  URL:  http://www.cecer.armv.mil. 
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2  Resuspension  of  Particles 

Resuspension  is  the  process  by  which  particulate  matter  is  introduced 
back  into  the  air  space  from  deposition  surfaces;  it  is  the  net  result  of  ad¬ 
hesion  forces  and  removal  forces  acting  on  particles  simultaneously  over 
time.  Resuspension  plays  a  significant  role  in  many  environmental  and  in¬ 
dustrial  processes  including:  transport  of  radionuclides,  exhaust  emission 
and  transport,  spreading  of  crop  diseases  by  fungal  spores,  cleaning  of 
electronic  chips,  handling  of  toxic  powders,  transmission  of  human  dis¬ 
eases,  and  distribution  of  allergens  in  indoor  spaces  (Gomes  et  al.  2005; 
Biasi  et  al.  2001). 

Early  scientific  studies  of  resuspension  processes  were  done  in  the  field  of 
erosion  and  soil  transport.  Development  of  nuclear  technology  led  to  re¬ 
search  on  resuspension  of  radionuclides  as  a  cause  of  spreading  airborne 
radioactivity  from  nuclear  weapon  test  sites  or  from  possible  accidental 
release  from  the  nuclear  industry.  Other  studies  in  resuspension  include 
the  transport  of  hazardous  particles  from  industrial  spills  and  pesticides 
(Nicholson  1988).  Recent  interest  in  resuspension  of  particles  in  indoor 
environments  has  been  raised  as  an  important  factor  in  indoor  air  quality 
(Thatcher  and  Layton  1995;  Hu  et  al.  2005). 

The  anthrax  attacks  on  U.S.  Postal  Service  facilities  and  the  Hart  Senate 
Office  Building  in  2001  demonstrated  that  the  release  of  a  small  amount  of 
CB  agents  in  a  building  is  a  very  effective  method  of  terrorism  that  poses  a 
serious  threat  to  the  building  occupants’  safety  and  that  can  disable  critical 
governmental  functions.  Airborne  CB  agents  released  in  one  section  of  a 
building  travel  via  the  building’s  heating,  ventilating,  and  air-conditioning 
(HVAC)  systems  and  disperse  throughout  the  building,  while  undergoing 
deposition  and  resuspension  recurrently. 

To  provide  proper  protection  from  CB  agents,  to  design  countermeasure 
devices,  and  to  plan  decontamination  schemes,  it  is  important  to  under¬ 
stand  how  CB  agents  migrate  through  a  building.  CFD  has  been  employed 
to  simulate  dispersion  of  CB  agents  in  a  building  through  HVAC  systems. 
How  well  the  simulation  result  would  match  the  actual  dispersion  profile 
of  CB  agents  depends  on  the  accuracies  of  models  describing  the  physical 
phenomena  such  as  aerodynamic  characteristics  of  CB  agents  and  fluid 
motions  through  HVAC  systems.  One  of  the  important  factors,  not  well  es- 


ERDC/CERL  TR-08-4 


5 


tablished,  is  the  resuspension  mechanism  of  fine  particles.  This  is  due  to 
the  complexity  of  the  turbulent  boundary  layer  where  particles  reside  and 
the  chemistry  of  particle  and  surface  interactions  (Loosmore  2003).  Re¬ 
suspension  models  currently  used  are  simple  and  inadequate  for  describ¬ 
ing  short-term  resuspension,  which  is  important  in  indoor  environments. 
Therefore,  a  resuspension  model  that  characterizes  indoor  resuspension 
appropriately  is  needed  that  would  improve  the  accuracy  of  simulation  of 
CB  agent  dispersal  in  a  building. 

2.1  Resuspension  Mechanism 

Resuspension  is  a  net  result  of  adhesion  forces  and  removal  forces  exerting 
concurrently  on  particles  deposited  on  a  surface.  This  section  examines 
adhesion  and  removal  forces. 

2.1.1  Adhesion  Forces 

The  main  adhesion  forces  between  a  particle  and  a  surface  are  the  van  der 
Waals  force,  the  capillary  force,  and  the  electrostatic  force. 

2. 1.1.1  The  van  der  Waals  Force 

Equation  1  describes  (and  Figure  2  shows)  the  van  der  Waals  force  be¬ 
tween  a  particle  and  a  near  surface  without  considering  contact  angle  flat¬ 
tening: 


F  = 

1  vdW 


^132^p 

12s2 


Eq  1 


where: 

A132=  the  Hamaker  constant 
dp  =  the  particle  diameter 

s  =  the  separation  distance  between  the  particle  and  the  surface. 


Figure  2.  Van  der  Waals  force  between  particle  and  plane. 


ERDC/CERL  TR-08-4 


6 


The  Hamaker  constant,  A132,  expresses  dipole  interactions  between  two 
materials  (Particle  l  and  Surface  2)  in  a  medium  (3,  air)  macroscopically: 

^132  =  ^12  ^33  _  ^13  _  ^23  ^ 

where: 

4  =  ,/44  Ec<3 

Combining  Equations  2  and  3  results  in: 

^132  =  (4/^H  _  V^33  ~  4/^33  )  4 

Real  particles  are  not  perfectly  rigid,  but  deform  under  the  influence  of 
this  adhesive  force.  If  the  adhesion  surface  area  is  increased  from  point 
contact  to  a  circle  of  radius  V  adhesion  ?  then  the  van  der  Waals  adhesion  force 
increases  to: 


_  ^132^p  ^  ^^132 ^adhesion 

12s2  12s3 


Eq  5 


Equation  6  describes  (and  Figure  3  shows)  van  der  Waals  force  between 
two  planar  surfaces: 


F  A 

1  vdW  _  ^132 

Area  6s3 


Eq  6 


Surface  1 

~  J"j  Medium  (Air)  3 

Surface  2 

Figure  3.  Van  der  Waals  force  between  planes. 

2. 1.1.2  Capillary  Force 

Depending  on  the  relative  humidity  (RH)  level,  water  can  condense  in  the 
gap  between  a  particle  and  a  surface  resulting  in  an  adhesion  force  caused 
by  surface  tension.  The  capillary  force  is  given  as: 

Fc  =  27ryLVdp  Eq  7 


where: 

yix  =  the  liquid-gas  surface  tension. 
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2. 1.1.3  Electrostatic  Force 

The  electrostatic  force  due  to  electrostatic  double  layer  is  expressed  as: 


77  —  Trf  — — 

1  EDL 


dP  <t>c 


Eq  8 


where: 

</>c  =  the  contact  potential 
so  =  the  dielectric  constant. 


In  case  of  charged  particles,  the  electrostatic  force  due  to  a  charged  parti¬ 
cle  is: 


^0  = 


1  <72 

4 7rs0  ( dp  +2s )2 


Eq  9 


where: 

q  =  the  charge  on  the  particle. 

2. 1.1.4  Empirical  Correlation 

Hinds  (1999)  derived  an  empirical  equation  for  adhesion  force  for  hard 
particles  and  clean  surfaces: 

Faih=  I5dp  (0.5  +  0.0045  (%RH)) 


The  units  for  Fadh  and  dp  are  N  and  pm,  respectively.  The  RH  term  in 
Equation  10  reflects  the  importance  of  the  capillary  force  in  adhesion 
forces. 

2.1.2  Removal  Forces 

A  removal  force  is  required  to  dislodge  particles  from  a  surface.  A  removal 
force  can  be  mechanical,  centrifugal,  vibrational,  or  hydrodynamic.  This 
section  considers  only  the  hydrodynamic  force.  Particles  on  a  surface  are 
subjected  to  shear  stress  in  the  viscous  sub-layer  of  the  turbulent  bound¬ 
ary  layer,  which  results  in  the  lift  and  the  drag  forces. 
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2.1.2.1  Lift  Force 

The  shear  lift  originates  from  inertia  effects  in  the  viscous  flow  around  the 
particle.  The  expression  for  the  inertia  shear  lift,  first  obtained  by  Saffman 
(1965;  1968)  can  be  restated  for  a  particle  resting  on  a  surface  as: 


F,  ( du, 

—lt  =  0.807 
oy  \  v 


where: 

pa  =  the  density  of  air 
v  =  the  kinematic  viscosity  of  air 
u  =  the  friction  velocity. 

Alternatively,  Leighton  and  Acrivos  (1985)  found: 


F,  ( dnii, 

— F-  =  0.576 
oy  \  v 


Eq  11 


Eq  12 


Hall  (1988)  derived  an  empirical  correlation  for  the  mean  lift  force  as: 


hr  *20.9 
2  2v 


2.1.2.2  Drag  Force 


d  u  ^ 

for  1.8  <  — e —  <  70 
2v 


Eq  13 


The  drag  force  for  a  particle  on  a  surface  is  given  by  Punjrath  and 
Heldman  (1972): 


C,  p  U2 1  nd] 

fx  r  i\  00  P 


Eq  14 


where: 

Cfx  =  the  local  shear  stress  coefficient 
Uoo  =  the  free  stream  air  flow  velocity. 

Hall  (1988)  derived  an  empirical  correlation  for  the  mean  drag  force: 


Eq  15 
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2.2  Existing  Models  of  Resuspension 

2.2.1  Resuspension  Factor  K  and  Resuspension  Rate  A 

The  most  frequently  used  terms  for  describing  resuspension  are  the  resus¬ 
pension  factor  K  and  the  resuspension  rate  A.  The  resuspension  factor  is 
defined  as: 

„  r  n  airborne  concentration  [g  m"3 1 

A:[m  ]  = - - - - — 

initial  surface  concentration  [g  m'  ]  Fn  If 


The  airborne  concentration  is  measured  at  some  reference  height,  which  is 
usually  the  breathing  height.  The  resuspension  rate  A  is  defined  as  the 
fraction  removed  in  unit  time,  more  explicitly: 

^  !  _  resuspension  flux  [g  m'V1  ] 

initial  surface  concentration  [g  m"2  ]  Fn  1 7 


The  relationship  between  K  to  A  is  given,  assuming  the  particle  concentra¬ 
tion  follows  a  power  law  distribution  with  height  by  (Loosmore  2003): 


A  =  ufhpK 


with: 


U»K 


where: 


Eq  18 


Eq  19 


k  =  the  von  Karman  constant 

Va  =  the  turbulent  deposition  velocity 

The  diabatic  influence  function  <p  has  a  value  of  1  for  neutral  conditions,  <1  for  unsta¬ 
ble  conditions,  and  >1  for  stable  conditions. 

The  value  for  p  has  been  calculated  to  be  in  the  range  of  0.05  to  0.5  with  a  median 
value  of  0.18. 


2.2.2  Empirical  Models 

Table  1  lists  empirical  models  for  the  resuspension  factor  K,  which  de¬ 
creases  exponentially  or  inversely  with  time. 
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Table  l.  Models  for  the  resuspension  factor  K. 


Model* 

if  =  1 0  5  exp(-0.0 1 390  + 1 0^9 
(Source:  USAEC  1972) 

Eq  20 

K  =  1 0"5  exp(-0.00 1 85r)  + 1  O'9 

Eq  21 

(Source:  USAEC  1975) 

K  =  10“4  exp(-0.15y7)  +  10”9 

Eq  22 

(Anspaugh  et  al.  1975) 

i:  =  10^exp(-0.010  +  10'9 

Eq  23 

(Linsley  1978) 

^  1.2  xlO-6 

K  = 

t 

(Garland  1979) 

Eq  24 

1  rr6 

K  =  10'6  for  t  <\  day  ;  K=  for  \<t  <1000  days 

(NCRP  1999) 

Eq  25 

K  =  [1 0^5  exp(-0.070  +  6  x  1 0~9  exp(-0.0030  + 1 0~9  ]  x  1 0±[ 

Eq  26 

(Anspaugh  et  al.  2002) 

K  =  — —  for  t  <  1000  days 

(Loosmore’s  (2003)  adaptation  of  NCRP  1999) 

Eq  27 

*  The  units  for  K  are  m"1  and  t  is  in  days. 

Loosmore  tested  the  applicability  of  the  existing  models  (Equations  22,  25, 
26,  and  27  in  Table  1)  to  short-term  resuspension  occurring  within  the  first 
24  hrs  and  found  that  only  Equation  27  is  adequate  to  describe  short-term 
resuspension  (Loosmore  2003).  The  author  also  developed  two  empirical 
models,  one  extensive  and  the  other  simplified,  for  short-term  resuspen¬ 
sion  using  published  data.  The  extensive  model  (Loosmore  2003)  is  ex¬ 
pressed  as: 


2.13  7  0.17 

u* 

A  =  0.42-  p 


0.92  0.32  0.76 


"0 


Po 


where: 

/7p  =  the  particle  density 

z0  =  the  aerodynamic  surface  roughness  length. 


Eq  28 
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The  simplified  model  developed  by  Loosmore  is  shown  as: 


where  (for  Equations  28  and  29): 
A  is  in  s-1 
u*  is  in  m/s 
t  is  in  s 
dp  is  in  pm 
Pp  is  in  kg/m3 
z0  is  in  m. 


Eq  29 


Lengweiler  et  al.  concluded  that  the  air  velocity,  turbulence  level,  and  sur¬ 
face  orientation  play  a  major  role  in  deposition  and  resuspension  (Leng¬ 
weiler  et  al.  1998).  Lengweiler  derived  empirical  correlations  address  that 
the  resuspension  rate  A  [s1]  is  proportional  to  turbulent  kinetic  energy  k 
[m2/s2]  for  floor,  wall,  and  ceiling  separately  (Lengweiler  2000).  They  are: 


A  =  1.85x10"4A: 

A  =  2.37xKr4£  +  1.08xl(T4 
A  =  3. 94x1 0“4  A: +  1.47x1 0”4 


(for  floor) 

Eq  30 

(for  wall) 

Eq  31 

(for  ceiling) 

Eq  32 

2.2.3  Theoretical  Models 


Wen  and  Kasper  (1989)  developed  a  kinetic  ‘particle  desorption’  model  in 
analogy  to  the  desorption  of  molecules  from  heterogeneous  surfaces.  Par¬ 
ticles  on  the  surface  have  a  distribution  of  adhesion  forces;  the  total  parti¬ 
cle  density  on  the  surface  is  obtained  by  integrating  the  particle  density  (as 
a  function  of  adhesion  force)  over  the  entire  range  of  the  adhesion  force: 

/•GO 

Nt  =  f  N(F)dF 

J»  Eq  33 

The  change  in  the  number  density  of  particles  with  the  adhesion  force  F  on 
the  surface  with  time  is  given  as: 


dN(F) 

dt 


=  - aN(F ) 


where: 


Eq  34 


N(F)  =  the  number  density  of  particles  with  the  adhesion  force  F  on 
the  surface. 

a  =  the  rate  constant. 
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The  first  order  rate  constant  is  a  function  of  the  non-dimensional  adhesion 
force,  F,  and  analogous  to  the  Arrhenius  formula: 


a(F)  =  Ae~F 

with: 

^  adhesion  force 

F  = - 

removal  force 


Eq  35 
Eq  36 


For  F,  Wen  and  Kasper  (1989)  used  a  uniform  initial  distribution.  Jurcik 
and  Wang  (1991)  developed  instead  a  kinetic  model  with  the  concept  of 
distributions  of  adhesion  and  removal  forces.  They  assumed  that  adhesion 
forces  have  a  log-normal  distribution  and  removal  forces  have  a  Gaussian 
distribution. 

The  above-mentioned  kinetic  models  are  also  called  force  balance  models. 
These  models  are  assuming  that  there  is  a  threshold  for  removal  based  on 
a  balance  between  the  instantaneous  aerodynamic  removal  force  and  par¬ 
ticle  surface  adhesive  forces.  Once  the  threshold  is  exceeded,  the  rate  of 
removal  is  determined  by  the  frequency  of  the  turbulent  bursts. 

Braaten  developed  a  force  balance  model  based  on  a  Monte-Carlo  ap¬ 
proach  instead  of  a  kinetic  approach  (Braaten  1994).  His  model  assumes 
that  aerodynamic  forces  are  applied  to  the  particle  on  the  surface  in  dis¬ 
crete  bursts.  The  values  of  these  forces  are  derived  from  some  prescribed 
probability  distribution;  resuspension  occurs  when  the  force  exceeded  the 
local  adhesive  force,  which  was  given  a  log-normal  distribution.  This  kind 
of  model  is  heavy  on  computing  time,  but  offers  a  direct  way  of  imple¬ 
menting  the  physical  process  via  the  particle  equation  of  motion  (Biasi  et 
al.  2001). 

Reeks  et  al.  (1988)  described  resuspension  as  a  process  whereby  particles 
reside  in  potential  wells  formed  by  attractive  and  repulsive  forces.  The  par¬ 
ticles  are  also  subject  to  lift  forces  with  mean  and  fluctuating  components 
generated  by  turbulent  flow  over  the  surface.  Resuspension  occurs  when 
the  particle  acquires  enough  energy  to  escape  from  the  potential  well; 
therefore,  this  model  is  also  called  a  dynamic  energy  accumulation  model. 
Recently,  Reeks  and  Hall  (2001)  developed  the  dynamic  “rock  and  roll” 
model  with  resonant  energy  transfer.  This  model  accounts  for  influence  of 
drag  through  accumulation  of  rotational  energy  through  the  rocking  of  a 
particle  about  the  surface  asperities  in  the  contact  zone. 
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2.3  Responsible  Parameters  for  Resuspension 

Resuspension  is  governed  by  adhesion  and  removal  forces.  Parameters 
that  affect  these  forces  include:  the  size  and  the  shape  of  particle,  the  ma¬ 
terial  characteristics  of  particles,  the  roughness  of  the  surface,  the  material 
characteristics  of  the  surface,  the  friction  velocity,  the  turbulent  level  of  the 
flow,  the  relative  humidity,  the  temperature,  and  the  duration  of  contact 
(Hinds  1999). 

As  shown  in  Section  2.1  (p  5),  there  is  an  agreement  between  different 
models,  whether  theoretical  or  empirical,  that  in  general,  adhesion  forces 
are  proportional  to  dp  while  removal  forces  are  proportional  to  dp2  or  dp 3. 
This  indicates  that,  as  the  size  of  particles  decreases,  adhesion  forces  be¬ 
come  greater  than  removal  forces.  Corn  and  Stein  (1965)  showed  experi¬ 
mentally  that  removal  of  particles  increased  with  increasing  particle  size. 
The  dependence  of  removal  efficiency  on  particle  size  was  partly  explained 
by  the  smallest  particles  becoming  enclosed  in  the  viscid  boundary  later 
and  not  being  easily  removed  (Nicholson  1988).  Larger  particles  stick  up 
higher  into  the  boundary  layer  and  are  subject  to  higher  removal  forces 
(Loosmore  2003). 

Corn  and  Stein  (1965)  reported  that  particle  adhesion  increased  with  in¬ 
creasing  relative  humidity,  probably  due  to  increasing  capillary  force.  The 
relative  humidity  also  has  an  effect  on  the  electrostatic  force.  Higher  rela¬ 
tive  humidity  results  in  faster  leakage  of  electrostatic  charge  thus  changing 
the  charge  on  particles. 

The  equations  in  Section  2.1  for  adhesion  and  removal  forces  pertain  to 
smooth  surfaces;  they  do  not  take  the  surface  roughness  into  account.  Ad¬ 
hesion  forces  decrease  with  increasing  surface  roughness,  which  is  due  to 
increasing  separation  distance.  However,  experimental  data  showed  that 
the  removal  of  particles  decreased  with  increasing  surface  roughness 
(Corn  and  Stein  1965).  Loosemore  (2003)  explained  that  larger  surface 
roughness  provides  more  shielding,  which  acts  against  resuspension. 

Time  dependence  of  resuspension  has  been  reported  in  many  studies.  The 
resuspension  flux  is  shown  to  decrease  nonlinearly  with  time  after  deposi¬ 
tion.  Nicholson’s  (1993)  wind  tunnel  experiments  showed  that  nearly  50 
percent  of  the  total  resuspension  occurred  within  the  first  10  seconds  in 
many  trials.  This  is  believed  to  be  due  to  the  fact  that  particles  with  the 
smallest  adhesive  forces  are  thought  to  be  removed  immediately  after 
deposition,  leaving  behind  a  more  firmly  held  contaminant  reservoir 
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(Clever  and  Yates  1973;  Ziskind  et  al.  1995;  Zimon  1969).  Linsley  described 
that  the  resuspension  factor  decreases  exponentially  with  time  (Linsley 

1978) .  Garland  found  that  an  inverse  power  law  fit  his  wind  tunnel  data  for 
times  varying  from  several  minutes  to  a  number  of  months  (Garland 

1979) .  It  is  widely  accepted  that  the  resuspension  factor  is  proportional  to 
i/fn  where  n  is  some  power  slightly  greater  than  1. 
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3  Development  of  Indoor  Resuspension 
Models  of  Fine  Particles 

3.1  Identification  of  Relevant  Parameters 

The  following  parameters  were  identified  as  important  for  the  resuspen¬ 
sion  rate  A  [s_1]  based  on  the  responsible  mechanisms  for  particle  resus- 


pension: 

dP: 

particle  diameter  (m) 

/¥ 

particle  density  (kg/m3) 

t : 

time  (s) 

u< 

friction  velocity  (m/s) 

turbulent  intensity  (turbulent  energy  dissipation  rate)  (m2/s3) 

Pt: 

turbulent  pressure  (kg/m/s2) 

z0: 

surface  roughness  length  (m) 

A132: 

Hamaker  constant  (kg  m2/s2) 

RH: 

relative  humidity 

/LV- 

liquid-gas  surface  tension  (kg/s2) 

Pa- 

air  density  (kg/m3) 

Pa: 

air  viscosity  (kg/m/s) 

9- 

gravity  (m/s2). 

Particle  properties  are  the  diameter  and  the  density;  the  surface  roughness 
length  characterizes  the  surface.  Friction  velocity,  turbulent  intensity,  and 
relative  humidity  are  environmental  conditions.  There  should  also  be  a  pa¬ 
rameter  or  parameters  that  differentiate  types  of  particles  (e.g.,  spores, 
Ti02,  silica,  etc.)  and  types  of  surfaces  (e.g.,  carpet,  wood,  concrete,  vinyl 
etc.).  The  Hamaker  constant  is  regarded  here  to  be  one. 

This  study  also  considered  thermophoresis,  which  describes  the  motion  of 
suspended  particles  induced  by  temperature  gradients  in  fluid.  Thermo¬ 
phoresis  may  play  a  role  in  dispersion  of  particles  (and  CB  agents)  in  a 
building  through  HVAC  systems,  especially  in  summer  or  winter  when  air 
conditioning  or  heating,  respectively,  is  used.  It  may  also  enhance  or  re¬ 
duce  particle  deposition  if  the  temperature  difference  between  walls  and 
the  ambient  is  considerable.  However,  the  effect  of  the  difference  between 
the  ambient  temperature  and  the  wall’s  temperature  on  deposited  particles 
and  the  magnitude  of  the  resulting  forces  compared  to  van  der  Waals 
forces  are  unknown.  Therefore,  neither  the  temperature  nor  the  tempera¬ 
ture  difference  in  the  above  was  included. 
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3.2  Dimensional  Analysis 

3.2.1  Case  1 

With  all  the  variables  listed  above,  the  resuspension  rate  is  expressed  as  a 
function  of  13  variables: 

A  =  f(dp,pp,t,  u»,  s,  RH,  z0,  Am,  pa,  Pt,  g,/iy)  Eq  37 


In  Equation  37,  there  are  a  total  of  14  parameters  and  three  repeating 
units:  mass,  length,  and  time.  Using  the  Buckingham  n  theorem,  the 
minimum  number  of  unique  dimensionless  parameters  is  11;  dp,  u*,  and  pA 
were  chosen  as  the  three  repeating  parameters  that  encompass  all  three 
dimensions,  and  the  resulting  equation  becomes: 


AJp 


u * 


/ 


Pp_ 

P a 


«•'  £dV  RTT  Z0 

~T’  — r>  KH’  ~j~ ’ 
d„  u,  dn 


dvu,p,  ’ 


Pt  g^P  7lv 

2  ’  2  5  r  2 

p a  dvu*  p a 


Eq  38 


3.2.2  Case  2 


To  derive  empirical  correlations  from  Equation  38  requires  experimental 
data  that  contain  information  on  all  10  independent  variables.  However, 
most  resuspension  experiments  in  literature  do  not  specify  values  of  s  and 
pt.  Therefore,  it  is  not  possible  to  include  s  andpt  as  variables  for  resus¬ 
pension  rate. 


Some  of  variables  in  Equation  37  do  not  change  or  are  relatively  constant. 
The  liquid-gas  surface  tension  of  interest  is  only  the  surface  tension  of  wa¬ 
ter  at  the  air  interface  for  the  particle  resuspension.  Therefore,  the  value  of 
Yly  does  not  vary  among  experiments,  and  can  be  excluded  in  Equation  37. 
Since  the  goal  of  this  work  is  to  develop  an  indoor  resuspension  model,  the 
range  of  ambient  temperature  is  relatively  small.  Furthermore,  most  of  ex¬ 
periments  were  performed  between  20  and  25  °C.  Since  the  viscosity  of  air 
is  primarily  a  function  of  temperature  with  a  minor  dependency  on  pres¬ 
sure,  the  variation  of  the  air  viscosity  would  be  very  small  and  can  be  omit¬ 
ted  from  Equation  37. 


Reeks  and  Hall  (2001)  noted  that,  in  practice,  the  contribution  from  grav¬ 
ity  is  much  smaller  than  that  from  the  drag  force,  although  it  may  be  im¬ 
portant  in  re-establishing  contact  with  the  surface  once  the  contacts  have 
been  broken  (Ziskind  2006).  Therefore,  Equation  37  can  disregard  gravity. 
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After  considering  these  aspects,  Equation  37  can  be  rewritten  as: 

A  =  f(d  p  ,pp,t,u*,KH,z0,AU2,pa) 


Eq  39 


and  the  resulting  dimensionless  variables  by  choosing  dp,  u*,  and  pa  as  the 
three  repeating  parameters  are: 


Adp 


u » 


f 


V 


P  a 


A 


Eq  40 


3.2.3  Case  3 

The  density  of  air  is  a  function  of  pressure  and  temperature.  In  indoor  en¬ 
vironments,  the  ambient  pressure  and  the  temperature  are  relatively  con¬ 
stant.  Therefore,  the  above  reasoning  for  considering  the  viscosity  of  air 
constant  can  be  applied  to  the  density  of  air,  and  the  relevant  parameters 
can  be  further  reduced  for  the  resuspension  as: 


A  =  f(dp ,pp,t,u » ,  RH,  Zq ,  Am ) 


Eq  41 


with: 


dp,  u*,  and  /y>,  instead  of  pa,  as  the  three  repeating  parameters,  thus: 


Ad 


-  =  / 


132 


Utt  DU  Z0 

“T’  ’  tt-t 

v  p 


dp  dp  u*  p 


Eq  42 


py 


3.2.4  Case  4 

To  be  thorough,  this  case  was  added  to  include  the  viscosity  of  air  to  Equa¬ 
tion  39.  As  a  result: 


A  =  /(dp,  Pp,t,  u* ,  RH,  z0 ,  Al32 ,  pa,  pa) 


Using  dp,  u*,  and  pA  as  the  three  repeating  parameters  again: 

f  •  1 


A  d„ 


=  f 


P  P  U*t  UTT  Z0 

- ,  —r~ ,  Kn,  — , 

P„  d„  d„ 


A 


132 


Pi 


dp  u*  P&  dpu*pa 


The  last  dimensionless  variable  in  Equation  44  is  l/Re. 


Eq  43 


Eq  44 


3.3  Model  Development 

Two  different  routes  were  pursued  for  model  development:  (1)  to  develop 
empirical  correlations,  and  (2)  to  derive  a  model  based  on  physics  involved 
in  resuspension  using  a  moment  balance  between  adhesion  and  removal. 
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3.3.1  Empirical  Correlations 

3.3. 1.1  Experimental  Data  Used  To  Develop  Empirical  Correlations 

Experimental  data  were  collected  from  references  by  Ibrahim  (2004), 
Nicholson  (1993),  Giess  et  al.  (1997),  Braaten  et  al.  (1990),  Wu  et  al. 

(1992),  Krauter  and  Biermann  (2007),  and  Wen  and  Kasper  (1989). 
Among  these  data,  Ibrahim’s  were  the  most  extensive  and  encompassed 
effects  of  the  relative  humidity  and  different  particle  sizes,  densities,  and 
materials.  However,  only  one  surface  material  (glass)  was  used  for  study¬ 
ing  resuspension.  Therefore,  another  data  set  from  a  study  using  different 
surface  materials  was  needed.  Consequently,  this  work  used  the  data  by 
Ibrahim  (2004)  and  Nicholson  (1993)  to  obtain  empirical  correlations. 

The  data  were  not  used  for  deriving  empirical  correlations,  but  were  used 
to  test  the  performance  of  the  empirical  correlations. 

Table  2  lists  the  parameters  used  for  the  experimental  works  by  Ibrahim 
and  Nicholson.  One  distinctive  characteristic  of  Ibrahim’s  experiments 
was  that  the  air  flow  was  transient,  meaning  the  velocity  of  the  flow  line¬ 
arly  accelerated  from  o  m/s  to  a  set  value,  whereas  most  of  experiments  in 
literature  used  a  constant  free  steam  velocity.  In  addition,  Ibrahim  also 
used  larger  (dp  =  ~  70  pm)  and  heavier  particles  (pp  =  8000  kg/m3)  than 
those  used  in  most  of  literature. 

Ibrahim  measured  the  surface  height  distribution  of  glass  using  atomic 
force  microscopy  (AFM).  The  standard  deviation  of  the  height  ranged  from 
1.08  nm  to  13.75  nm  with  the  mean  of  the  standard  deviation  of  1.7  nm. 
Because  the  physical  roughness  or  the  effective  obstacle  height  is  usually 
estimated  by  the  difference  between  the  lowest  and  the  highest  points  for  a 
surface,  the  highest  value,  13.75  nm,  was  taken  as  the  standard  deviation 
value  for  the  glass  surface.  If  one  can  assume  the  height  has  a  normal  dis¬ 
tribution,  99.7  percent  of  the  values  are  within  three  standard  deviations. 
Therefore,  the  physical  roughness  was  calculated  as  three  times  the  stan¬ 
dard  deviation.  The  physical  roughness  of  grass  in  Nicholson’s  data  was 
given.  Loosmore  (2003)  assumed  the  physical  roughness  for  concrete  in 
Nicholson’s  to  be  3  mm,  and  this  work  used  the  same  value.  Loosmore  es¬ 
timated  the  aerodynamic  roughness  length  as  1/10  of  the  physical  rough¬ 
ness.  Krauter  and  Biermann  (2007)  used  1/30  of  the  physical  roughness  as 
the  aerodynamic  roughness  based  on  the  work  by  Sutton  (1955).  Here, 

1/10  of  the  physical  roughness  was  used  as  the  aerodynamic  roughness 
length. 
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Table  2.  Experimental  parameters  for  works  by  Ibrahim  (2004)  and 

Nicholson  (1993). 


Ibrahim  (2004) 

Nicholson  (1993) 

Dimensions  for  flow  passage 

Wind  tunnel,  0.203  m  x  0.203  m 

Wind  tunnel,  lmxlm 

Surface  material 

Glass 

Grass  and  concrete 

Uco  (m/s) 

Transient,  0-20 

For  grass:  3.0,  4.5,  8.0 

For  concrete:  3.0,  5.0,  6.5 

u*  (m/s) 

u*  =  0.03751/co  +  0.0387 

For  grass:  0.15,  0.3,  0.4 

For  concrete:  0.1,  0.15,  0.2 

Particle 

Stainless  steel  (SS)  spheres, 
glass  spheres,  Lycopodium 
spores 

Silica  spheres 

dp  (urn) 

SS:  70,  glass:  72,  32,  Lyco¬ 
podium:  30 

4.1,  9.6,  17.5,  22.1  for  most 
exp.’s 

4.8,  8.6,  12.1,  16.5  for  1  set 

Pe  (kg/m3) 

SS:  8000,  Glass:  2420,  Lyco¬ 
podium:  1000 

1000 

Surface  feature  height  (m) 

4.1x10-8 

Grass:  0.05,  concrete:  0.003 

Aerodynamic  roughness  (m) 

4.1x10-9 

Grass:  0.005,  concrete: 
0.0003 

Temperature  (°C) 

23 

Not  reported,  assumed  to  be 
25 

RH  (%) 

Varied, 18  -  67 

Not  reported 

Resuspension  rate 

Computed  from  detachment  frac¬ 
tion  using  Figs.  5.2,  5.3,  5.4,  5.7, 
5.9,  6.2,  6.5,  8.1  in  Ibrahim 
(2004) 

Taken  from  Fig.  2  in  Nichol¬ 
son  (1993) 

Number  of  data  points 

443 

96 

Ibrahim’s  data  (2004)  were  presented  as  the  detachment  fraction  vs.  the 
free  stream  velocity.  The  velocity  linearly  accelerated  during  the  experi¬ 
ments  as  mentioned  earlier,  and  therefore  the  velocity  was  converted  to 
time  using  the  acceleration  value.  Once  graphs  were  obtained  as  the  de¬ 
tachment  fraction  vs.  time,  then  resuspension  rate  was  calculated  by  tak¬ 
ing  a  slope  between  two  adjacent  points.  For  the  corresponding  time  for 
the  obtained  resuspension  rate,  two  different  approaches  can  be  used  as 
discussed  by  Nicholson  (1993)  and  Loosmore  (2003).  The  simplest  ap¬ 
proximation  for  the  corresponding  time  is  to  use  the  mid-point  of  the  pe¬ 
riod.  An  alternative  approach  is  to  assign  the  representative  time  to  reflect 
the  fact  that  the  resuspension  rate  falls  inversely  proportionally  to  the 
time.  The  corresponding  time  t'  for  the  calculated  resuspension  rate  can  be 
obtained  by: 


Eq  45 
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Here,  it  was  decided  to  use  both  the  mid-point  and  t  to  draw  empirical 
correlations  and  compare  the  results.  Nicholson’s  data  were  taken  directly 
from  Figure  2  in  Nicholson  (1993),  in  which  the  mid-point  time  was  used 
for  the  corresponding  time  for  the  measured  resuspension  rate.  Since  the 
sampling  intervals  were  given  in  the  reference,  t'  was  also  calculated  as  for 
Ibrahim’s  data.  That  is,  these  two  data  sets  encompass  the  time  range  of 
1  -  10000  s  and  the  resuspension  rate  range  of  io-8  -  101  s_1  (Figure  4). 


Time  (s) 

Figure  4.  Resuspension  rate  versus  time  plotted  from  the  data  of  Ibrahim 
(2004)  and  Nicholson  (1993).  For  the  corresponding  time,  (a)  used  the  mid¬ 
point  and  (b)  used  t'  defined  by  Eq  45. 
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The  Hamaker  constants  for  Ibrahim’s  and  Nicholson’s  data  were  calcu¬ 
lated  using  Equation  4.  For  example,  to  calculate  the  Hamaker  constant 
for  stainless  steel  (particle)-glass  (surface)-air  (medium),  the  Hamaker 
constants  for  stainless  steel-stainless  steel-air  (vacuum)  and  glass-glass-air 
(vacuum)  are  needed. 

Table  3  lists  the  Hamaker  constants  for  homogenous  substances  that  were 
used  for  Equation  4.  The  Hamaker  constant  for  stainless  steel-stainless 
steel-vacuum  was  assumed  to  be  the  same  as  that  of  steel-steel-vacuum. 
The  Hamaker  constant  for  Lycopodium- Lycopodium-air  was  calculated 
using  Equation  4  and  the  Hamaker  constants  for  Lactobacillus-glass- 
water,  water-water-air,  and  glass-glass-air,  assuming  that  biological  mate¬ 
rials  have  similar  values  for  the  Hamaker  constant  since  the  molecular 
compositions  are  similar.  Under  the  same  assumption,  the  value  of  the 
Hamaker  constant  for  grass-grass-vacuum  was  considered  as  the  same  as 
that  of  Lycopodium-Lycopodium-air  once  obtained.  The  Hamaker  con¬ 
stant  for  concrete-concrete-air  is  not  available,  and  it  was  assumed  to  be 
5x10-2°  J. 

3.3. 1.2  Effect  of  RH 

As  mentioned  earlier,  Ibrahim’s  data  were  the  most  extensive  with  varia¬ 
tions  in  relative  humidity,  particle  sizes,  densities,  and  materials.  Table  4 
briefly  lists  experimental  data  from  Ibrahim  (2004);  and  SS70  was  used  in 
the  most  of  experiments.  If  all  the  results  using  SS70  are  pooled  (Figures 
5.2,  5.3,  5.4,  5.7,  5.9,  and  8.1  from  Ibrahim  [2004]),  the  particle  size,  the 
density,  and  materials  for  particle  and  surface  are  constant  for  these  data 
sets.  Variables  for  these  experiments  are  RH,  the  acceleration,  and  the  fi¬ 
nal  free  stream  velocity. 


Table  3.  Hamaker  constants  used  for  calculating  A132  for  Ibrahim’s  and 

Nicholson’s  data. 


Material 

A  (J) 

Reference 

steel-steel-vacuum 

21.2x10-20 

Soltani  and  Ahmadi  (2004) 

glass-glass-dry  air 

8.5x10-20 

Soltani  and  Ahmadi  (2004) 

SiC>2-Si02-vacuum 

6.82x10-2° 

French  (2000) 

Lycopodium-Lycopodium-air 

6.5x10-20 

Calculated  from  Lactobacillus-glass-water  (6.2x10- 
21  J,  Rijnaarts  et  al.  (1995)),  water- water-air 
(3.7x10-2°  J),  and  glass-glass-air 

grass-grass  -air 

6.5x10-20 

Assumed  to  be  the  same  as  Lycopodium- 
Lycopodium-air 

cone  rete-co  n  c  rete-a  i  r 

5x10-20 

Assumed 
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Table  4.  Variables  for  the  experimental  work  by  Ibrahim  (2004). 


Figure 

No. 

from 

Ibrahim 

2004 

Particle 

Particle 

Density 

(kg/m3) 

Surface 

T 

(°C) 

Acceleration 

(m/s2) 

Initial 

Particle 

No. 

Density 

(mm-2) 

RH 

(%) 

5.2 

SS70 

8000 

Glass 

23 

0.18 

0.5 

Varied 

Effect  of  RH 

5.3 

SS70 

8000 

Glass 

23 

Varied 

0.5 

36 

Effect  of  acceleration 

5.4 

SS70 

8000 

Glass 

23 

Varied 

0.5 

61 

Effect  of  acceleration 

5.7 

SS70 

8000 

Glass 

23 

0.18 

0.5 

30 

Effect  of  final  free  stream 
velocity 

5.9 

SS70 

8000 

Glass 

23 

0.15 

0.5 

25 

Variability  of  results  for 
the  same  experimental 
condition 

6.2 

GL72 

GL32 

2420 

Glass 

23 

0.18 

0.5 

25 

Effect  of  particle  size 

6.5 

GL32 

LY30 

2420 

1000 

Glass 

23 

0.15 

0.5 

25 

Effect  of  particle  density 

8.1 

SS70 

1000 

Glass 

23 

0.014 

0.5 

52 

Variability  of  results  for 
the  same  experimental 
condition 

The  final  free  stream  velocity  did  not  affect  the  fraction  of  particle  resus¬ 
pended,  as  noted  by  the  author,  and  can  be  eliminated  from  the  variables. 
In  contrast,  the  flow  acceleration  during  the  transient  period  did  have  an 
effect  on  the  fraction  of  particle  resuspended.  However,  the  acceleration  is 
not  the  parameter  considered  here  for  resuspension,  and  it  was  decided  to 
neglect  the  effect  of  the  flow  acceleration.  Consequently,  the  pooled  data 
for  SS70  has  only  one  variable  (RH).  Note  that  the  data  that  were  dupli¬ 
cated  among  Ibrahim’s  figures  were  deleted  while  being  pooled. 

For  the  pooled  data  with  SS70,  the  particle  size,  the  particle  density,  the 
surface  roughness  length,  and  the  Hamaker  constant  were  constant.  Then, 
Equation  39  can  be  simplified  further  to  for  these  pooled  data  as: 

A  =  /(f,n*,  RH,  pa)  Eq  46 


And,  the  resulting  dimensionless  variables  are: 


M, 

w* 


f 


=f\ 


V 


u*  t 


\ 

,  RH 
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Eq  47 
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Assuming  a  power  law  for  each  variable,  Equation  47  becomes: 


A  d 


=  a  0 


c  v1 

Uih  t 


d 

p  J 


(l  +  a2RHfl3) 


Eq  48 


It  can  also  be  expressed  as: 


fA0 

(  \ 

£/*  t 

In 

=  ln(aO)  +  a\  •  In 

J 

cl 

\  V  J 

+  ln(l  +  a2-RH"3) 


Eq  49 


Equation  49  implies  that  the  plot  of  \n(Adp/u*)  vs.  \n(u*t/dp)  at  a  specific 
RH  is  a  straight  line  and  the  slope  is  at.  In  addition,  the  slope  does  not 
change  at  different  RH,  but  the  intercept  does. 

Figure  5  shows  the  graph  of  \n(Adp/u*)  vs.  \n(u*t/dp)  at  various  RH.  Most 
of  the  data  form  a  cluster  around  what  could  be  a  straight  line  with  some 
scatter  regardless  of  RH.  Moreover,  the  scatter  is  not  confined  to  a  certain 
RH,  but  is  noticeable  from  all  RH.  This  graph  does  not  illustrate  that  dif¬ 
ferent  intercepts  are  needed  for  different  values  of  RH,  and  does  indicate 
that  the  effect  of  RH  is  rather  small  and  can  be  neglected  in  the  range  of  18 
-  67  percent.  Therefore,  it  was  decided  to  exclude  the  RH  from  the  rele¬ 
vant  parameters  in  Equations  39,  41,  and  43. 


Figure  5.  Plot  of  hif/lc/p/ii  )  vs.  ln(u*t/dp)  of  all  SS70  results  with  different  RH 

from  Ibrahim  (2004). 
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3.3. 1.3  Empirical  Correlations 

If  one  assumes  that  a  power  law  applies  for  each  dimensionless  number  in 
Equation  40  and  the  effect  of  RH  is  negligible,  then  a  general  correlation 
can  be  expressed  as: 


M, 
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=  b0 


Pp_ 

\P*J 


vM 


,b  2 


U *  t 
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(  -  \ 
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Eq  50 


Taking  the  logarithms  of  both  sides  yields: 
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Eq  51 


To  obtain  the  coefficients  60-65,  a  multiple  regression  was  done  using 
Ibrahim’s  and  Nicholson’s  data  with  the  mid-point  as  the  corresponding 
time  for  the  resuspension  rate.  Table  5  lists  the  results  of  a  multiple  re¬ 
gression,  and  the  correlation  becomes: 
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Eq  52 


Table  5.  The  multiple  regression  result  of  Eq  51  using  Ibrahim’s  (2004)  and 
Nicholson’s  (1993)  data  with  the  mid-point  as  the  corresponding  time  for  the 

resuspension  rate. 


R  =  0.96046136  R2= 

0.92248603  Adjusted  R2  =  0.92190540 

Standard  Error  of  Estimate  =  1.0451 

Coefficient 

Standard  Error 

t 

P 

In(foO) 

-4.7652 

0.7346 

-6.4865 

<0.0001 

bl 

-0.3028 

0.0855 

-3.5397 

0.0004 

bl 

-1.0135 

0.0263 

-38.5816 

<0.0001 

b3 

-0.3269 

0.0200 

-16.3183 

<0.0001 

M 

-0.2961 

0.0383 

-7.7326 

<0.0001 

Analysis  of  Variance: 

DF 

ss 

MS 

F 

p 

Regression 

4 

6941.3353 

1735.3338 

1588.7701 

<0.0001 

Residual 

534 

583.2612.3711 

1.0922 

Total 

538 

7524.5967 

13.9862 

ERDC/CERL  TR-08-4 


25 


The  R2  value  for  the  correlation  is  0.92,  and  the  probabilities  for  any  coef¬ 
ficient  to  be  zero  is  very  low  (<0.0001  -  0.0004).  The  exponent  for  the  t 
term  is  -1.0145,  consistent  with  the  findings  of  Nicholson  (1993)  and 
Loosmore  (2003). 

The  same  procedures  were  used  to  derive  empirical  correlations,  denoted 
EC  I  to  EC  VII,  for  different  groups  of  dimensionless  variables  in  Cases  2 
through  4  (Sections  3.2.2-3.24,  pp  16-17)  while  assuming  the  effect  of 
RH  on  the  resuspension  rate  was  negligible.  Each  dimensionless  function 
was  evaluated  using  both  the  mid-point  and  t  for  the  corresponding  time 
for  the  resuspension  rate;  Table  6  lists  the  results.  Although  all  the  statis¬ 
tics  of  each  multiple  regression  are  not  shown,  the  probability  of  the  F-test 
was  <0.0001  for  all  empirical  correlations.  The  probabilities  for  any  coeffi¬ 
cient  to  be  zero  from  the  two-tailed  t-test  were  also  very  low;  in  most  of 
cases  p<o.oooi  and  in  all  cases  p<0.05  except  for  EC  VI.  In  EC  VI,  p  from 
the  t-test  for  ln(do)  and  cLj.  was  0.2387  and  0.0713,  respectively.  There¬ 
fore,  EC  VI  was  re-evaluated  while  letting  ln(do)  =  o  and  d4  =  o;  the  re¬ 
sulting  correlation  is  shown  as  EC  VII  (which  does  not  contain  the 
Hamaker  constant  as  a  variable). 

The  data  listed  in  Table  6  show  that,  for  the  same  group  of  dimensionless 
numbers,  choosing  either  the  mid-point  or  t'  as  the  representative  time  for 
the  resuspension  did  not  affect  the  results  of  empirical  correlations  signifi¬ 
cantly;  the  exponents  for  each  dimensionless  parameter  are  similar  and  R2 
values  are  changed  little.  A  better  fit  was  expected  using  t'  as  the  corre¬ 
sponding  time  for  the  resuspension,  but  the  mid-point  time  showed  very 
slightly  better  fits  for  all  three  different  groups  of  dimensionless  numbers. 

The  empirical  correlations  (EC  I  -  EC  VII)  were  non-dimensionalized  to 
examine  the  dependency  of  A  on  the  variables;  Table  7  lists  the  results.  For 
all  empirical  correlations,  the  exponent  of  the  time  term  is  between  -0.92 
and  -1.04;  these  values  are  similar  to  those  in  Loosmore’s  work  (2003)  or 
the  findings  by  Nicholson  (1993).  All  these  correlations  suggest  that  A  in¬ 
creases  with  u*  to  the  power  n  with  n<i,  and  this  is  smaller  than  that  sug¬ 
gested  by  Loosmore.  Sehmel  also  noted  that  the  air  concentration  is  found 
to  increase  with  wind  speed  to  power  n,  with  n  between  1  and  6.42  (Ri- 
jnaarts  et  al.  1995).  Therefore,  the  exponents  for  u*  in  the  correlations  here 
seem  too  low  compared  to  values  reported  in  the  literature.  On  the  other 
hand,  the  dependency  on  the  particle  diameter  in  these  correlations  is 
stronger  than  that  of  Loosmore  (2003). 
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Table  6.  Empirical  correlations  obtained  using  different  groups  of 
dimensionless  numbers.  Each  dimensionless  function  was  evaluated  using 
both  the  mid-point  and  t'  for  the  time  that  represents  the  resuspension  rate, 
and  the  effect  of  the  RH  was  assumed  to  be  negligible. 


with  mid-point 


with  t' 


From  Equation  40 


Ad 


=  b0 


\P*J 


UA*'-  Y3 


Y4 


\d*y 


\dpy 


132 


Kdv  u*  P*J 


EC  I 

bO  =  8.521x10-3 
bl  =-0.3028 
b2  =-1.0135 
b3  =-0.3269 
b4  =-0.2961 
R2  =  0.9224 


EC  IV 

bO  =3.753x10-3 
bl  =-0.2921 
b2  =-0.9712 
b3  =-0.3646 
b4  =-0.2820 
R2  =  0.9215 


From  Equation  42 


Ad 


—  =  c0 


v  c3 


\d*j 


\dvj 


1132 


i  3  2 

dp  u,  Pv 


ecu 

c0  =  3.708x10-4 
cl  =-0.9635 
c2  =-0.3224 
c3  =-0.1742 
R2  =  0.9172 


EC  V 

c0  =1.882x10-4 
cl  =-0.9255 
c2  =-0.3576 
c3  =-0.1672 
R2  =  0.9158 


From  Equation  44 


Ad 


-  =  d0 


/ 

P  ' 


\P»  j 


f  y1 

u*  t 


vrfpy 


( .  Y1 


4 


,d  5 


EC  III 

C/0  =  5.193x102 
dl  =-0.4096 
d2  =-1.0959 
d3  =-0.3408 
d4  =  0.4815 
d5  =-1.8636 
R2  =  0.9249 


Y32 


3  2 

d p  u *  p a 


ydpU*pa  j 


EC  VI 

c/0  =  2.509x101 
dl  =-0.3761 
d2  =-1.0336 
d3  =-0.3791 
d4  =  0.3468 
d5  =-1.500 
R2  =  0.9231 
EC  VII 
dO  =  1 
dl  =-0.4644 
d2  =-1.0352 
d3  =-0.3649 
d4  =  0 
d5  =-0.7861 
R2  =  0.9221 


Table  7.  Dependency  of  A  on  each  variable  for  different  models. 


Model 

Exponent  of 

t 

u* 

cfp 

Pp 

zo 

A132 

//a 

EC  1 

-1.0135 

0.5787 

1.2311 

-0.3028 

-0.3269 

-0.2961 

0.5989 

- 

ecu 

-0.9635 

0.3849 

0.8085 

0.1742 

-0.3224 

-0.1742 

- 

- 

EC  III 

-1.0959 

0.8047 

0.8558 

-0.4096 

-0.3408 

0.4815 

1.7917 

-1.8636 

EC  IV 

-0.9712 

0.5928 

1.1818 

-0.2921 

-0.3646 

-0.2820 

0.5741 

- 

EC  V 

-0.9255 

0.4089 

0.7847 

0.1672 

-0.3576 

-0.1672 

- 

- 

EC  VI 

-1.0336 

0.7728 

0.8723 

-0.3761 

-0.3791 

0.3468 

1.5293 

-1.5000 

EC  VII 

-1.0352 

0.7509 

1.1862 

-0.4644 

-0.3649 

- 

1.2505 

-0.7861 

Eq  28* 

-0.92 

2.13 

0.17 

-0.76 

-0.32 

- 

- 

- 

Eq  29  * 

-1.03 

1.43 

- 

- 

- 

- 

- 

- 

*  from  Loosmore  (2003) 
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EC  I  and  EC  IV  illustrate  that  A  increases  with  increasing  u*  and  dp  and 
decreases  with  increasing  t,  pp,  and  z0.  This  is  consistent  with  the  findings 
by  others,  and  the  reasoning  for  these  dependencies  is  described  in  Section 
2-3  (p  13)  and  can  also  be  found  elsewhere  (Loosmore  2003).  In  addition, 
EC  I  and  EC  IV  also  implied  that  A  increases  with  decreasing  A132  and  in¬ 
creasing  pa.  The  lower  value  of  the  Hamaker  constant  implies  that  the  van 
der  Waals  forces  between  particles  and  surfaces  are  smaller,  and  thus  the 
particles  can  be  easily  resuspended.  The  higher  value  of  air  density  would 
increase  the  forces  from  the  flow  to  the  particles  and  increase  the  resus¬ 
pension  rate. 

EC  II  and  EC  V  were  the  correlations  using  the  dimensionless  numbers  in 
Equation  40,  where  pa  was  eliminated  from  the  relevant  parameters.  These 
two  correlations  indicate  that  the  resuspension  rate  increases  with  increas¬ 
ing  particle  density.  Loosmore  noted  that  the  gravitational  attraction  is 
expected  to  be  much  smaller  than  the  van  der  Waals  or  capillary  adhesion 
forces  and,  therefore,  the  particle  density  is  not  believed  to  affect  the  re¬ 
suspension  rate  directly  (Loosmore  2003).  However,  the  author  also  men¬ 
tioned  that  deposition  or  re-deposition  after  resuspension  depends 
strongly  on  particle  mass  and,  thus,  particle  density.  Therefore,  particles 
with  a  higher  density  would  result  in  a  lower  net  resuspension  rate  for  the 
same  size  particles,  which  is  the  opposite  of  predictions  by  EC  II  and  EC  V. 
Since  they  render  physically  unrealistic  models,  EC  II  and  EC  V  are  not 
considered  further. 

EC  III  and  EC  VI  include  pa  as  a  parameter  and  suggest  that  the  resuspen¬ 
sion  rate  increases  with  decreasing  the  fluid  viscosity  or  increasing  Re.  Al¬ 
though  this  is  reasonable,  the  models  also  indicate  that  the  resuspension 
rate  increases  as  the  Hamaker  constant  increases.  A  higher  Hamaker  con¬ 
stant  means  higher  van  der  Waals  adhesion  and  should  result  in  a  lower 
resuspension  rate.  Therefore,  EC  III  and  EC  VI  do  not  make  physical  sense 
and  are  not  considered  further. 

In  summary,  only  EC  I,  EC  IV,  and  EC  VII  (which  does  not  include  the 
Hamaker  constant  as  a  variable)  are  acceptable. 

3.3.2  Physics-Based  Model 

A  physics-based  model  was  also  developed  using  a  moment  balance  be¬ 
tween  adhesion  and  removal.  According  to  literature,  the  most  realistic 
way  of  particle  resuspension  is  by  rolling  or  tangential  lift.  Therefore,  the 
condition  for  particle  detachment  was  assumed  to  be  when  the  removal 
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moment  is  bigger  than  the  adhesion  moment.  In  another  words,  if  the 
moment  ratio,  M,  is  defined  as: 


Adhesion  Moment 
Removal  Moment 


Eq  53 


then  M<1  is  the  resuspension  condition.  For  the  adhesion  force,  using  the 
Johnson-Kandall-Roberts  (JKR)  model  (Johnson  et  al.  1971),  the  adhesion 
moment  is  expressed  as: 


Ma  =  — ;zrcr-r 

u  2  P  3 


Eq  54 


Here,  rp  is  the  particle  radius,  cr  is  the  surface  energy,  and  ra  is  the  contact 
radius.  The  hydrodynamic  removal  moment  is  from  the  drag  force,  and 
therefore: 


Mr=\AFDrv 


Eq  55 


FD  =  6mjrnf 


Eq  56 


Here,  Y  is  the  shear  rate  and/=  1.7009.  Because  the  friction  velocity  and 
the  shear  rate  are  related  as: 


f  x 

/  •  \ 
Aa  r 

UaJ 

V  Pa  J 

Eq  57 


Equation  53  becomes: 

,,  1  arA 

M  = - 5 — 

5.6/  w.2pa  r2 


Eq  58 


and  Equation  58  can  be  rearranged  with  two  dimensionless  numbers  as: 

..if,  Yr.’l 


V 

r  \ 

a 

2 

u*  p  r 

V  #P  y 

r 

v  p  y 

Eq  59 
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or  with  thee  dimensionless  numbers  as: 


M  =■ 


5.6/ 


Y.  ^ 


u*  Pp  rP 


r 

v  p; 


V^a  J 


Eq  60 


The  contact  radius,  ra,  at  the  moment  of  separation  is  given  by: 


r  = 


'aO 


a  /l  1/3 


Eq  61 


where  ra0  is  the  contact  radius  at  zero  applied  load  given  as: 


ra0  = 


^pvv/3 

V  *  , 


Eq  62 


and 


3 


1-v,  l-v22 

1  +-  2 


Ex  E2 


-i 


Eq  63 


Here,  £i  and  E2  are  Young’s  moduli  for  the  particle  and  the  surface,  respec¬ 
tively,  and  Vi  and  v2  are  Poisson  ratios  of  the  particle  and  the  surface,  re¬ 
spectively.  The  resuspension  rate,  A,  is  defined  as: 


A  =  - 


1  clN 
N0  dt 


Eq  64 


where: 

No  and  N  are  the  number  density  of  particles  on  the  surface  at  time  o 
and  t,  respectively. 

Wen  and  Kasper  (1989)  assumed  a  Langmuir  model  for  the  particle  de¬ 
sorption  kinetics.  Using  the  same  assumption: 

N  =  N0e~at  Eq  65 

and  the  rate  constant  a  is  a  function  of  M. 

According  to  Wen  and  Kasper,  a  requirement  for  a(M)  is  to  decrease  with 
increasing  M.  The  two  simplest  functions  that  meet  the  requirement  are: 

a  =  Axe~M 


Eq  66 
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M 


By  choosing  Equation  66,  Equation  64  can  be  rewritten  as: 

A  =  A{e~u-e~Ae~M‘ 

By  introducing  the  definition: 


and  rearranging  Equation  65,  one  obtains: 

A  =  —e~tlT' 

h 

Alternatively,  choosing  Equation  67  results  in: 

A  =  —e~tlTl 
T 2 

with: 

M 


Eq  67 


Eq  68 


Eq  69 


Eq  70 


Eq  71 


Eq  72 


For  given  particles,  M  is  a  constant.  Therefore,  a  model  can  be  obtained  by 
fitting  Equation  70  or  71  to  experimental  data  to  find  A  or  A2,  respectively. 

3.3.3  Comparison  of  Dimensionless  Variables  in  Empirical  Correlations 
and  Physics-based  Model 

One  of  the  motivations  for  developing  a  physics-based  model  was  to  de¬ 
termine  dimensionless  variables  that  are  based  on  physics  and  examine 
whether  they  can  be  used  for  drawing  empirical  correlations.  The  dimen¬ 
sionless  numbers  from  the  physics-based  model  in  Equation  59  were  com¬ 
pared  to  those  used  for  empirical  correlations  in  Equation  40.  In  Equa¬ 
tion  59,  two  dimensionless  numbers  take  account  for  the  properties  of 
particles  and  surfaces:  <j/u*2parv  and  ra/rp.  In  Equation  40,  the  term  with 
the  Hamaker  constant  (Ai32/ dp3u*2pa)  takes  into  account  the  material 
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characteristics  of  particles  and  surfaces.  If  the  “contact”  interaction  from 
the  JKR  model  is  used  instead  of  the  van  der  Waals  interaction  for  the  em¬ 
pirical  correlations,  then  Equation  40  can  be  replaced  with: 


A  d 


=  f 


£  DU  ^0_  17 

Pa  dp  dp  u*  pad  d 


Eq  73 


p  J 


Here,  da  is  the  contact  diameter,  2 ra.  Acquiring  an  empirical  correlation 
with  the  dimensionless  numbers  in  Equation  73  will  be  pursued  in  Phase 
II  of  this  project. 

3.4  Testing  of  Models 

Empirical  correlations,  EC  I,  EC  IV,  and  EC  VII,  were  compared  against 
the  experimental  results  by  Giess  et  al.  (1977),  Braaten  et  al.  (1990),  Wu  et 
al.  (1992),  Krauter  and  Biermann  (2007),  and  Wen  and  Kasper  (1989). 

3.4.1  Collected  Experimental  Data 

Figure  6  shows  the  data  from  the  literature  to  test  empirical  correlations; 
Tables  8  and  9  list  their  experimental  parameters. 


The  friction  velocity  for  the  data  by  Wu  et  al.  (1992)  was  calculated  using 
the  following  equation  for  duct  flows  (suggested  by  Davies  [1972])  and  the 
free  stream  velocity: 

5 


c  V/8 

v 


ut  = 


V^h^oo  J 


Eq  74 


where: 


dh  =  the  hydraulic  diameter. 


Using  work  by  Krauter  and  Biermann  (2007)  and  Wen  and  Kasper  (1989) 
the  following  equation  for  a  pipe  flow  was  used  to  calculate  the  friction  ve¬ 
locity  (Davies  1966): 


where: 

f  =  the  friction  factor  given  by  Blasius  as: 


0.3164 
4  Re025 


Eq  76 


ERDC/CERL  TR-08-4 


32 


O 

'to 

c 

0) 

a 

to 

3 

to 

<D 

□£ 


03 

oc 


(/) 

C 

<u 

Q. 

CO 

3 

10 

<D 

0£ 


IO1 
10° 

10-1 
io-2 
io-3 
io4 

IQ-5 
IQ-8 
IQ 7 

10-1  100  101  10?  103  104  105 

Time  (s) 


— i — i  i  i  1 1  n| - 1 — . | - 1 — i  . . . . 1 — »  i  i  rin| - 1 — i  i  imi| - l — . . 


A  g  A 

I 

&  o 


■  ' 


A 


o  Giess 
r  □  Braaten 
a  Wlj 
.  O  Krauter 
0  Wen 

_ ■  ■  ■  ■  ^  _ *  . . . ■  *  » ■■'■J _ i _ l 


A 

A  A 


O  ] 


101 
10° 

10-1 
io-2 
io-3 
10^ 

105 
io-6 
io-7 

10*1  10°  101  102  103  104  105 


a 

A 

A 

A 


A 
,  A 


V 

S3**- * 

° 


B) 


18 


Giess 

Braaten 

Wu 

Krauter 

Wen 

-mil. 


Time  (s) 

Figure  6.  Resuspension  rate  versus  time  plotted  from  the  data  used  to  test 
empirical  correlations.  For  the  corresponding  time,  (A)  used  the  mid-point 

and  (B)  used  t  defined  by  Eq  45. 


Table  8.  Experimental  parameters  for  works  by  Giess  (1997)  and 

Braaten  (1990). 


Giess  (1997) 

Braaten  (1990) 

Dimensions  for  flow  passage 

Wind  tunnel,  1.2  m  x  0.8  m 

Wind  tunnel,  lmxlm 

Surface  material 

Grass 

Glass 

ur_  (m/s) 

3.0, 5.0, 7.8 

6.0,  7.5,  9.0 

u*  (m/s) 

For  short  grass:  0.27,  0.64,  1.14 
For  long  grass:  0.57,  0.75,  1.42 

0.17,  0.215,  0.26 

Particle 

Silica  spheres 

Lycopodium  spores 

c/p  (|xm) 

1.85 

27.8 

Pp  (kg/m3) 

1000 

1000 

Surface  feature  height  (m) 

Short  grass:  0.1 

Long  grass:  0.3 

Not  reported 
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Giess  (1997) 

Braaten  (1990) 

Aerodynamic  roughness  (m) 

Short  grass:  0.01 

Long  grass:  0.03 

1x10-6;  assumed  by 
Loosmore  (2003) 

Temperature  (°C) 

Not  reported,  assumed  to  be  25 

Not  reported,  assumed 
to  be  25 

RH  (%) 

Not  reported 

Not  reported 

Resuspension  rate 

Taken  directly  from  Table  1 

Computed  from  fraction 
loss  of  experimental  data 
in  Fig.  12  -  Fig.  14 

Number  of  data  points 

30 

21 

Table  9.  Experimental  parameters  for  works  by  Wu  (1992),  Krauter  (2007), 

and  Wen  (1989). 


Wu  (1992) 

Krauter  (2007) 

Wen  (1989) 

Dimensions  for  flow 
passage 

Wind  tunnel,  lmxlm 

Round  duct,  D  =  0.152  m 

Capillary  tubing, 

D  =  0.125  cm 

Surface  material 

Glass 

Steel,  plastic 

Steel 

Uo o  (m/s) 

4,  6,8 

2.6 

100 

u*  (m/s) 

0.169,  0.241,  0.310: 
calculated 

0.144:  calculated 

5.16:  calculated 

Particle 

Lycopodium  spores 

Bacillus  atrophaeus 
spores 

Latex  particles 

dp  (pm) 

27.8 

0.91 

0.509,  1.019 

Pp  (kg/m3) 

1000 

1200 

1000 

Surface  feature 
height  (m) 

Not  reported 

Steel:  1.5X10-4, 
plastic:5xl0-6 

Not  reported 

Aerodynamic 
roughness  (m) 

lxlO'6:  assumed  by 
Loosmore6 

Steel:  5xl0  6, 
plastic:1.67xl0-7 

lxlO'6:  assumed 

Temperature  (°C) 

Not  reported,  assumed 
to  be  25 

Steel:  27,  averaged 

Plastic:  24,  averaged 

22 

RH  (%) 

Not  reported 

Steel:  26,  averaged 

Plastic:  26,  averaged 

~  0 

Resuspension  rate 

Computed  from  the 
resuspension  fraction 
from  Fig.  4 

Taken  directly  from  Fig.  6 

Computed  from  the 
particle  concentration 
in  air  in  Figs.  8  and  9 

Number  of  data  points 

96 

54 

81 

Other  details 

Pressure:  6  bar  p&: 
7.18  kg/m3 

The  aerodynamic  roughness  for  the  data  by  Giess  (1997)  is  assumed  to  be 
1/10  of  the  physical  roughness  as  discussed  earlier.  For  the  data  by 
Braaten  et  al.  (1990)  and  Wu  et  al.  (1992),  Loosmore  assumed  lxio-6  m  as 
the  aerodynamic  roughness  for  glass  surfaces;  it  was  decided  to  use  the 
same  value  for  this  study.  For  reference,  the  aerodynamic  roughness  for 
glass  surface  used  by  Ibrahim’s  work  was  estimated  as  4.1x10-9  m.  As  men- 
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tioned  earlier,  Krauter  and  Biermann  (2007)  estimated  the  aerodynamic 
roughness  as  1/30  of  the  physical  height  of  the  surface;  and  values  were 
used  directly.  The  surface  roughness  for  the  steel  capillary  tubing  used  by 
Wen  and  Kasper  (1989)  was  assumed  to  be  lxio-6  m. 

The  resuspension  rates  that  were  calculated  from  the  resuspension  frac¬ 
tion  used  the  mid-point  time  as  the  corresponding  time  for  EC  I  and  t'  for 
EC  IV  and  EC  VII  for  consistency.  The  resuspension  rate  of  Wen  and 
Kasper  was  computed  from  the  particle  concentration  in  air  using  the  flow 
rate,  the  initial  particle  number  density  on  the  surface,  and  the  surface 
area  deposited  with  particles.  Table  10  lists  the  homogenous  Hamaker 
constants,  in  addition  to  those  in  Table  3,  for  the  experimental  data. 

3.4.2  Model  Testing 

Figure  7  shows  the  ratio  of  the  resuspension  rate  predicted  by  empirical 
correlations  to  the  experimental  results.  The  performances  of  the  empiri¬ 
cal  correlations  of  EC  I,  EC  IV,  and  EC  VII  were  similar  to  all  experimental 
data.  For  the  data  by  Braaten  et  al.  (1990)  and  Wu  et  al.  (1992),  the  ratios 
of  the  predicted  to  the  observed  were  mostly  scattered  between  0.1  and  10 
for  all  three  correlations.  The  ratios  of  the  predicted  to  the  observed  by 
Wen  and  Kasper  (1989)  were  distributed  between  0.01  and  1  for  EC  I  and 
EC  IV.  EC  VII  performed  slightly  better  than  EC  I  or  EC  IV  for  Wen  and 
Kapser  (1989),  and  this  was  the  most  distinguishable  difference  in  com¬ 
parison  among  the  correlations.  The  ratio  was  relatively  constant  at  ~o.oi 
for  Giess  et  al.  (1997),  and  the  correlations  greatly  underestimated  the  data 
by  Krauter  and  Biermann  (2007). 


Table  10.  Additional  Hamaker  constants  used  for  calculating  A132  for  Giess  et 
al.  (1997),  Braaten  et  al.  (1990),  Wu  et  al.  (1992),  Krauter  and  Biermann 
(2007),  and  Wen  and  Kasper  (1989). 


Material 

A  (J) 

Reference 

B.  atrophaeus-B.  atrophaeus-air 

6.5x10-20 

Assumed  the  same  as  Lycopodium-Lycopodium-air 

latex  sphere-latex  sphere-air 

4x10-20 

Assumed 

plastic-plastic-air 

7.8x10-2° 

Assumed  the  same  as  PVC-PVC-vacuum 
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Figure  7.  Performance  of  the  empirical  correlations,  (A)  EC  I,  (B)  EC  IV,  and 
(C)  EC  VII,  against  experimental  data. 
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Note  that  Loosmore’s  empirical  correlations  predicted  the  data  by  Braaten 
et  al.  (1990)  and  Wu  et  al.  (1992)  reasonably  well  (Loosmore  2003).  The 
work  by  Wen  and  Kasper  was  done  in  a  capillary  tube  with  a  very  high  free 
stream  velocity  (100  m/s);  the  data  may  not  be  comparable  to  others  ob¬ 
tained  with  low  velocities.  Krauter  and  Biermann  (2007)  compared  their 
experimental  data  to  Loosmore’s  empirical  correlation,  and  prediction  by 
Loosmore’s  correlation  was  lower  by  two  orders  of  magnitude  than  the  ex¬ 
perimental  data.  They  also  compared  their  data  to  those  from  the  litera¬ 
ture  by  adjusting  for  the  friction  velocity  difference  using  the  exponent  for 
u*  in  Loosmore’s  correlation.  The  resuspension  rates  from  Braaten  et  al. 
and  Wu  et  al.  were  an  order  of  magnitude  lower  than  those  by  Krauter  and 
Biermann,  and  only  the  results  by  Wen  and  Kasper  were  comparable. 

Because  values  for  the  surface  roughness  were  not  available  for  most  of  the 
experimental  works,  they  were  estimated  to  evaluate  the  empirical  correla¬ 
tions.  This  can  introduce  or  increase  uncertainties  in  the  empirical  correla¬ 
tions.  The  correlations  exponents  for  z0  were  between  -0.32  and  -0.36.  If 
the  surface  roughness  would  be  1/ 1000  of  that  specified  in  Table  8  and 
Table  9,  it  would  cause  an  order  of  magnitude  increase  in  the  resuspension 
rate.  Therefore,  inaccurate  estimation  of  the  surface  roughness  would  not 
affect  the  prediction  of  the  resuspension  rate  too  greatly. 

The  quality  of  an  empirical  correlation  depends  on  the  data  used  to  draw 
the  correlations,  and  some  of  the  experimental  details  on  the  data  by  Ibra¬ 
him  (2004),  which  was  used  in  conjunction  with  the  data  by  Nicholson 
(1993),  deserve  discussion.  Specifically,  Ibrahim’s  data  were  obtained  us¬ 
ing  transient  velocities  compared  to  constant  velocities  used  by  most  of  the 
other  works.  In  addition,  some  of  Ibrahim’s  data  were  collected  using  very 
dense  particles  ( pp  =  8000  kg/m3).  Since  the  empirical  correlations  were 
obtained  using  the  dimensionless  variables,  some  differences  may  have 
been  scaled  properly  and  the  correlations  can  be  valid  for  wide  ranges  of 
parameters.  However,  it  may  be  possible  that  other  differences  may  be 
fundamental,  and  empirical  correlations  derived  from  these  data  may  not 
properly  describe  the  resuspension  data  under  different  conditions.  There¬ 
fore,  to  model  and  predict  the  indoor  resuspension  accurately  it  would  be 
essential  to  derive  models  based  on  experimental  data  with  the  conditions 
that  are  characteristic  to  indoor  resuspension. 
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Table  11  lists  the  results  of  simple  sensitivity  testing  for  the  empirical  cor¬ 
relations  of  EC  I,  EC  IV,  and  EC  VII.  A  factor-of-io  error  in  dv  would  result 
in  a  factor  of  15-17  change  in  the  resuspension  rate  for  all  correlations.  It 
is  not  likely  to  have  a  factor-of-io  error  in  the  particle  size,  and  the  contri¬ 
bution  from  this  error  would  be  very  small.  An  order  of  magnitude  error  in 
u *,  pp,  z0,  A 13 2,  and  pa  would  cause  only  a  factor  of  2  -  4  change  in  A  except 
for  pa  in  EC  VII,  and  the  contribution  from  the  errors  in  these  parameters 
would  be  limited. 

Table  11.  Sensitivity  analyses  of  the  empirical  correlations  of 
EC  I,  EC  IV,  and  EC  VII. 


Factor-of-10 
changes  in  the 
parameter  of 

Change  in  A 

ECI 

EC  IV 

EC  VII 

t 

0.097 

0.107 

0.092 

u* 

3.791 

3.916 

5.635 

dp 

17.026 

15.198 

15.353 

Pp 

0.498 

0.510 

0.343 

zo 

0.471 

0.432 

0.432 

A 132 

0.506 

0.522 

- 

Pa 

3.971 

3.750 

17.803 

Pa 

- 

- 

0.164 
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4  Source  Term  Model  for  CFD  Simulation 

Demonstration  simulations  were  performed  using  the  empirical  correla¬ 
tion  EC  I  as  a  source  term  in  a  boundary  condition  for  CFD  simulations. 
The  CFD  simulations  were  based  on  the  experiments  performed  by  Law¬ 
rence  Berkeley  National  Lab;  a  description  of  the  experiments  follows. 

4.1  Description  of  Particle  Resuspension  in  a  Bench-top  Test 
Chamber 

A  30  in.  (L)  x  30  in.  (W)  x  36  in.  (H)  chamber  was  used  for  particle  resus¬ 
pension  experiments.  The  chamber  contained  two,  3-blade  muffin  fans 
(blade  span  =  4.25  in.,  115  V,  14  W,  variable  speed)  mounted  10  in.  from 
the  floor  and  located  in  opposite  corners  of  the  chamber.  The  fans  blew 
towards  the  center  of  the  chamber,  and  the  center  of  each  fan  was  6  in. 
from  its  respective  corner  (Figure  8).  Make-up  air  at  20  L/min  was  intro¬ 
duced  to  the  chamber  from  a  0.25  in.  diameter  port  in  the  center  of  the 
right  wall  of  the  chamber,  which  was  located  18  in.  above  the  floor  and 
15  in.  from  the  front  wall  and  projected  0.25  in.  into  the  chamber.  An  aero¬ 
dynamic  particle  sizer  (APS,  Model  3321,  TSI,  Inc.,  St.  Paul,  MN)  drew 
5  L/min  of  air  from  a  probe  at  the  center  of  the  floor.  The  probe  extended 
5.25  in.  above  the  floor  and  had  a  0.75  in.  diameter.  In  addition,  10  L/min 
of  air  was  sampled  for  filter  measurement  through  a  port  24  in.  above  the 
chamber  floor  and  10  in.  from  the  front  wall.  It  extended  5.75  in.  into  the 
chamber  and  had  a  0.25-in.  diameter.  At  5  L/min,  air  leaked  along  the  pe¬ 
rimeters  of  26  in.  (W)  x  32  in.  (H)  of  a  door  on  the  front  wall.  Horizontal 
and  vertical  components  of  air  velocity  were  measured  using  a  hotwire 
anemometer  at  9  points  in  a  plane  10  cm  above  the  floor  specified  in  Fig¬ 
ure  8. 

During  the  particle  deposition  phase,  the  airborne  particle  concentration 
was  measured  with  the  APS,  and  the  data  collected  by  the  APS  were  used 
to  estimate  the  size  distribution  of  particles  deposited  on  the  floor  sur¬ 
faces.  During  the  particle  resuspension  phase,  the  airborne  concentration 
was  measured  with  both  the  APS  and  a  filter  sample.  The  resuspension 
phase  for  each  experiment  lasted  5  hrs,  1  hr  with  the  fans  off  and  1  hr  for 
each  of  four  fan  speeds  (45,  75, 110,  and  140).  The  fan  speeds  were  incre¬ 
mentally  increased  each  hour  over  the  course  of  the  experiment. 
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5  L/min  Air  Leak 

along  with  the  perimeter 
of  26”  x  32" 


10  L/min  Air  Out 

5.75"  projection  into  the  chamber 
0.25"  diameter 


3.94"  from  the  floor 


Figure  8.  A  schematic  diagram  of  the  chamber  for  particle  resuspension 
experiments;  all  dimensions  are  in  inches,  and  the  x-y  positions  of  the 
measurement  points  for  velocity. 


The  floor  of  the  chamber  was  completely  covered  with  either  moderately 
textured  linoleum  or  commercial  grade  olefin  carpet,  and  the  sodium  fluo¬ 
rescein  powder  was  used  as  the  particle. 

4.2  CFD  Simulations  of  Flow  Field  in  the  Test  Chamber 

4.2.1  Experimental  Results 

The  data  from  the  experiments  that  provided  information  on  the  flow  field 
were  the  flow  rate  of  the  fans  and  the  average  horizontal  and  vertical  ve¬ 
locities  at  the  plane  10  cm  above  the  floor;  Tables  12, 13,  and  14  list  these 
values. 
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Table  12.  Measured  flow  rates  of  the  variable  speed  fan. 


Fan  Speed 

Flow  Rate 
(L/min) 

45 

6 

75 

20 

110 

40 

140 

70 

Table  13.  The  average  of  the  minimum,  the  average,  and  the  maximum  air 
velocity  in  the  horizontal  direction  measured  at  the  nine  measurement 
points  in  the  plane  10  cm  above  the  floor  for  the  four  different  fan  speeds. 


Fan  Speed 

Average  of 
Minimums  (cm/s) 

Average  of 
Averages  (cm/s) 

Average  of 
Maximums  (cm/s) 

45 

1.17 

7.80 

16.83 

75 

17.89 

30.10 

46.30 

110 

24.26 

37.27 

54.00 

140 

21.61 

38.60 

57.71 

Table  14.  The  average  of  the  minimum,  the  average,  and  the  maximum  air 
velocity  in  the  vertical  direction  measured  at  the  nine  measurement  points  in 
the  plane  10  cm  above  the  floor  for  the  four  different  fan  speeds. 


Fan  Speed 

Average  of 
Minimums  (cm/s) 

Average  of 
Averages  (cm/s) 

Average  of 
Maximums  (cm/s) 

45 

5.85 

14.00 

27.15 

75 

30.57 

56.93 

94.54 

110 

28.40 

54.62 

98.81 

140 

24.59 

57.96 

103.33 

4.2.2  CFD  Simulations  of  Flow  Field 

CFD  simulations  in  this  project  were  carried  out  in  collaboration  with  CFD 
Research  Corporation  (Huntsville,  AL).  Based  on  the  chamber  geometry 
and  flow  conditions,  a  computational  mesh  was  created  using  CFD- 
Micromesh,  a  mesh  generation  software  developed  and  commercialized  by 
CFDRC.  Figure  9  shows  model  geometry  and  mesh  distribution  used  in  the 
simulations.  Localized  mesh  refinement  was  used  in  regions  where  sharp 
flow  gradients  are  expected,  e.g.,  inlets,  outlets,  and  regions  affected  by  fan 
flow  (Figure  9[B]).  Three  different  mesh  sizes  were  used  to  ensure  mesh- 
independent  results  of  flow  field.  A  total  of  45,113  grid-cells  were  used  in 
the  parametric  study. 

Once  the  mesh  was  generated,  flow  field  was  calculated  inside  the  chamber 
based  on  the  conditions  provided  using  CFD  code  CFD-ACE+,  originally 
developed  by  CFDRC  and  currently  licensed  by  ESI  Group  (San  Diego, 
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CA).  Figure  10  shows  the  Graphical  User  Interface  demonstrating  the 
range  of  physical  models  available  in  CFD-ACE+.  The  “modules”  that  were 
used  in  the  current  simulation  were  “flow,”  “turbulence,”  and  “user  scalar.” 
A  “fan  sub-model”  was  used  for  simulating  flow  and  pressure  head  of  the 
fan.  The  fans  were  assumed  to  have  an  outer  radius  of  5.2  cm  and  a  hub 
radius  of  2.975  cm. 


Figure  9.  (A)  3D  model  of  chamber  showing  the  general  direction  of  air  flow 
through  various  inlets  and  outlets.  (B)  unstructured  (binary-tree)  mesh 
implemented  inside  the  chamber  3D  model.  The  refined  mesh  layer  just 
below  the  make-up  air  vent  is  to  resolve  the  flow  modification  through  the 

two  fans. 
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Figure  10.  Snapshot  of  the  CFD-ACE+  graphical  user  interface  with  physical 

models  used  in  the  simulations. 

Parametric  simulations  of  flow  were  performed  for  various  fan  flow  rates. 
For  comparison,  Table  15  lists  the  measured  and  calculated  fan  flow  rates, 
and  the  flow  field  generated  by  CFD  simulations  that  could  match  the  fan 
flow  rates  from  the  experiments.  For  reference,  the  average  fan  flow  veloc¬ 
ity  was  calculated  based  on  the  measured  fan  flow  rate  and  the  radius  of 
the  fan  (either  for  the  outer  area  or  the  hub  area);  the  last  two  columns  in 
Table  15  lists  these  values.  Table  16  lists  the  average  horizontal  and  vertical 
velocities  of  the  nine  points  in  the  plane  that  was  10  cm  above  the  floor. 
The  measured  average  velocities  were  taken  directly  from  Tables  13  and 
14,  and  the  calculated  values  were  the  results  of  CFD  simulations  of  the 
flow  field. 

That  data  in  Table  16  show  that  the  average  calculated  velocities  are  sig¬ 
nificantly  lower  than  measured  values.  Even  when  the  flow  rate  was  in¬ 
creased  to  132.7  L/min,  which  was  much  higher  than  any  of  the  experi¬ 
mental  values,  the  average  horizontal  and  vertical  velocities  were  still 
considerably  lower  than  any  of  experimental  values.  Because  one  could 
simply  calculate  the  average  fan  flow  velocity  from  the  fan  flow  rate  and 
the  annulus  area  of  the  fan  as  shown  in  the  last  two  columns  in  Table  15, 
researchers  also  examined  the  average  horizontal  and  vertical  velocities  in 
the  plane  that  contained  the  fans  (Table  17).  The  average  horizontal  veloc¬ 
ity  in  the  plane  10  in.  above  the  floor  was  slightly  higher  than  the  average 
fan  flow  velocity  (in  Table  15)  for  6-40  L/min  and  lower  for  70  and 
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132.7  L/min,  but  overall  one  could  conclude  that  these  two  values  agreed 
acceptably.  Therefore,  there  seemed  to  be  inconsistency  between  the 
measured  flow  rate  and  the  measured  average  velocities.  It  was  decided  to 
use  the  measured  fan  flow  rate  for  the  basis  of  CFD  simulation  on  particle 
resuspension. 


Table  15.  Measured  and  calculated  fan  flow  rates;  the  average  fan  flow 
velocities  were  calculated  from  the  measured  fan  flow  rates  from  the 
experiments  using  the  outer  radius  and  the  hub  radius  of  the  fans. 


Measured  Fan 
Flow  Rate 
(L/min) 

CFD 

Calculated 
Flow  Rate  for 
Fan  #1 
(L/min) 

CFD 

Calculated 
Flow  Rate  for 
Fan  #2 
(L/min) 

CFD 

Calculated 
Average  Fan 
Flow  Rate 
(L/min) 

Average  Fan 
Flow  Velocity  - 
Based  on 
Outer  Area 
(cm/s) 

Average  Fan 
Flow  Velocity  - 
Based  on 
Annulus  Area 
(cm/s) 

6 

6.56 

6.40 

6.48 

1.18 

1.75 

20 

20.36 

19.48 

19.92 

3.92 

5.83 

40 

40.39 

36.76 

38.57 

7.85 

11.67 

70 

61.68 

55.16 

58.42 

25.51 

37.91 

- 

150.36 

115.08 

132.7 

33.31 

49.52 

Table  16.  Measured  and  calculated  average  velocities  of  the  nine 
measurement  points  in  the  plane  10  cm  above  the  floor  for  various  fan  flow 
rates;  the  average  horizontal  and  vertical  velocities  are  from  Tables  13  and  14. 


Calculated 
(Measured)  Fan 
Flow  Rate 
(L/min) 

Average 

Horizontal  Velocity 
-  Measured 
(cm/s) 

Average  Vertical 
Velocity  - 
Measured 
(cm/s) 

Average 

Horizontal  Velocity 
-  Calculated 
(cm/s) 

Average  Vertical 
Velocity  - 
Calculated 
(cm/s) 

6.48  (6) 

7.80 

14.00 

5.03 

3.61 

19.92  (20) 

30.10 

56.93 

4.50 

4.14 

38.57  (40) 

37.27 

54.62 

4.10 

3.76 

58.42  (70) 

38.60 

57.96 

4.60 

4.06 

132.7 

- 

- 

9.79 

7.41 

Table  17.  Measured  and  calculated  average  velocities  of  the  nine  points  in 
the  plane  that  the  fans  were  placed  (10  in.  above  the  floor)  for  various  fan 
flow  rates;  the  x-y  positions  of  the  nine  points  were  the  same  as  those  of  the 
nine  measurement  points  in  Thl.  16,  but  the  z-positions  for  these  points  were 
10-in  from  the  floor  instead  of  10  cm. 


Calculated 
(Specified)  Fan  Flow 
Rate  (L/min) 

Average  Horizontal 
Velocity-Calculated 
(cm/s) 

Average  Vertical 
Velocity-Calculated 
(cm/s) 

6.48  (6) 

2.26 

1.91 

19.92  (20) 

11.00 

2.41 

38.57  (40) 

14.33 

1.89 

58.42  (70) 

18.67 

2.65 

132.7 

30.90 

5.12 
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Figures  11-15  show  the  flow  field  and  friction  velocities  in  the  chamber  at  a 
given  fan  flow  rate.  The  flow  field  is  plotted  in  the  plane  of  the  two  fans  as 
well  as  in  the  plane  joining  the  fan  centers  (shown  on  the  left).  The  figures 
on  the  right  show  isocontours  of  friction  velocities  at  the  floor  of  the 
chamber.  As  expected,  as  the  fan  flow  rate  is  increased  the  velocities  (in¬ 
cluding  the  friction  velocities),  and  the  turbulence  levels  in  the  chamber 
increase  significantly. 
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Figure  11.  (A)  Velocity  vectors  colored  by  magnitude  of  velocity  in  the  plane 
through  the  fan  centers;  vectors  are  also  plotted  in  the  plane  along  the 
diagonal  through  bottom;  isosurfaces  of  velocity  magnitude  in  the  plane  of 
the  two  fans  (6  L/min);  (B)  isocontours  of  the  friction  velocity  at  the  floor 

(6  L/min). 
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Figure  12.  (A)  Velocity  vectors  colored  by  magnitude  of  velocity  in  the  plane 
through  the  fan  centers;  vectors  are  also  plotted  in  the  plane  along  the 
diagonal  through  bottom;  isosurfaces  of  velocity  magnitude  in  the  plane  of 
the  two  fans  (20  L/min);  (B)  isocontours  of  the  friction  velocity  at  the  floor 

(20  L/min). 
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Figure  13.  (A)  Velocity  vectors  colored  by  magnitude  of  velocity  in  the  plane 
through  the  fan  centers;  vectors  are  also  plotted  in  the  plane  along  the 
diagonal  through  bottom;  isosurfaces  of  velocity  magnitude  in  the  plane  of 
the  two  fans  (40  L/min);  (B)  isocontours  of  the  friction  velocity  at  the  floor 

(40  L/min). 
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Figure  14.  (A)  Velocity  vectors  colored  by  magnitude  of  velocity  in  the  plane 
through  the  fan  centers;  vectors  are  also  plotted  in  the  plane  along  the 
diagonal  through  bottom;  isosurfaces  of  velocity  magnitude  in  the  plane  of 
the  two  fans  (60  L/min);  (B)  isocontours  of  the  friction  velocity  at  the  floor 

(60  L/min). 
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Figure  15.  (A)  Velocity  vectors  colored  by  magnitude  of  velocity  in  the  plane 
through  the  fan  centers;  vectors  are  also  plotted  in  the  plane  along  the 
diagonal  through  bottom;  isosurfaces  of  velocity  magnitude  in  the  plane  of 
the  two  fans  (132  L/min);  (B)  isocontours  of  the  friction  velocity  at  the  floor 

(132  L/min). 
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4.3  CFD  Simulation  of  Particle  Resuspension  and  Transport 

4.3.1  Resuspension  Rate 


As  mentioned  earlier,  the  empirical  correlation  EC  I  used  as  a  source  term 
in  a  boundary  condition  for  CFD  simulations  is  written  as: 


Adn  , 

P  _8  SOI  vin-3 

-0.3028 

(  \ 

U„t  1 

-1.0135 

(ecI 

-0.3269 

r  A  ) 

^132 

-0.2961 

—  O.JZ1  a  1VJ 

u* 

vft  J 

A) 

lrfp  j 

PI  3  2 

d  M*  0 

\  p  *  J 

Eq  77 


and,  therefore,  the  resuspension  rate  is  expressed  as: 
A  =  8.521xl0“3 


0.5787  ,  1.2311  0.5989 

u «  d  pa 


Eq  78 


P p 


132 


In  the  experiments,  the  floor  of  the  chamber  was  completely  covered  with 
either  moderately  textured  linoleum  or  commercial  grade  short  pile  olefin 
carpet.  Linoleum  is  made  from  solidified  linseed  oil  (linoxyn)  in  combina¬ 
tion  with  wood  flour  or  cork  dust  over  a  burlap  or  canvas  backing.  A  com¬ 
mercial  grade  olefin  carpet  is  made  of  polypropylene,  which  is  static  resis¬ 
tant  unlike  nylon  carpets.  The  physical  roughness  for  linoleum  and  the 
short  pile  carpet  were  assumed  to  be  6.4  pm  (Hallas  and  Shaw  2006)  and 
1  mm,  respectively,  and  the  aerodynamic  surface  roughness  was  taken  as 
1/10  of  the  physical  roughness.  The  particle  size  distribution  of  the  sodium 
fluorescein  powder  was  measured  using  the  APS,  and  it  ranged  between 
0.5  -  15  pm.  The  density  of  the  sodium  fluorescein  was  1.35  g/cm3,  and  the 
density  of  air  was  1.196x10-3  g/cm3  for  22  °C.  Table  18  lists  the  homoge¬ 
nous  Hamaker  constants  for  the  linoleum,  polypropylene,  and  sodium 
fluorescein.  The  assumed  value  of  the  Hamaker  constant  for  linoleum  was 
based  on  the  fact  that  typical  values  of  the  Hamaker  constant  for  con¬ 
densed  phases  (solid  or  liquid),  interacting  across  vacuum,  are  »io-^  J. 
When  hydrocarbons  are  treated  as  an  assembly  of  CH2  groups,  the  value 
obtained  for  it  is  around  sxicr20  J  (Martins  et  al.  2006).  The  Hamaker 
constant  for  sodium  fluorescein  was  also  assumed  to  be  sxio-20  J,  since  the 
molecules  are  composed  of  mostly  C,  H,  and  some  O. 


Table  18.  The  homogenous  Hamaker  constants  of  the  floor  materials  and  the 
particle  used  for  calculating  A132  for  the  chamber  experiments. 


Material 

^(j) 

Reference 

linoleum-linoleum-air 

5x10-20 

Assumed 

polypropylene-polypropylene-air 

7.1x10-20 

Martins  et  al.  (2006) 

Sodium  fluorescein-sodium  fluorescein-air 

7.8x10-20 

Assumed 
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4.3.2  CFD  Simulations  of  Particle  Transport  and  Resuspension 

After  the  flow  field  was  calculated,  particles  were  introduced  into  the  do¬ 
main  from  the  bottom  floor  as  a  function  of  the  friction  velocity.  The  parti¬ 
cle  concentration  was  calculated  assuming  that  particles  were  convected 
along  with  the  flow,  and  the  particle  diffusivities  were  calculated  using  the 
Stokes-Einstein  equation.  Table  19  lists  properties  of  particles  as  a  func¬ 
tion  of  the  particle  size.  The  diffusion  coefficient  is  small  for  particles  lar¬ 
ger  than  0.5  pm  as  compared  to  nano-sized  particles.  Gravitational  settling 
could  be  significant  for  large-sized  particles  as  seen  from  the  gravitation 
terminal  velocity.  However,  the  effect  of  gravity  was  neglected  in  the  simu¬ 
lations  during  Phase  I,  because  scalar  transport  formulation  cannot  ac¬ 
commodate  the  effect  of  body  forces,  such  as  gravity,  in  the  computational 
calculations.  One  can  employ  Lagrangian  particle  tracking  (“spray”  mod¬ 
ule)  approach  or  Eulerian-Eulerian  two-fluid  (“two-fluid”  module)  model 
to  include  all  the  relevant  forces  affecting  particle  transport.  Both  of  these 
capabilities  are  available  in  CFD-ACE+  and  will  be  employed  in  Phase  II 
work. 


Table  19.  Properties  of  particles  as  a  function  of  the  particle  size. 


Particle  Size 
(pm) 

Diffusivity 

(m2/s) 

Slip  Correction  Cc 

Terminal  Velocity 
(cm/s) 

Relaxation  Time 
(s) 

0.01 

5.39x10-8 

22.225 

9.08X10-8 

9.26xl0-9 

0.1 

6.95xl010 

2.867 

1.17xl0-4 

1.19xl0'7 

0.5 

6.45X1011 

1.330 

1.36xl0-3 

1.39xl0-6 

2.12 

1.23X1011 

1.077 

1.98xl0-2 

2.02x10-6 

4.04 

6.24xl012 

1.041 

6.94x10-2 

7.08x10-5 

10.3 

2.39xl012 

1.016 

4.41X101 

4.49xl0-4 

Particle  transport  subsequent  to  resuspension  was  calculated  for  the  fan 
flow  rates  of  20  and  132  L/min.  A  source  term  for  particle  was  calculated 
using  EC  I,  assuming  the  particle  size  was  0.5  pm  and  the  initial  deposi¬ 
tion  mass  of  the  particle  was  352  ng  (0.06  ng/cm2).  Linoleum  was  the  floor 
material.  The  resuspension  at  time  o  was  assumed  zero  and  resuspension 
started  at  1  second;  Figure  16  shows  the  particle  source  from  resuspension. 
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Time  (s) 

Figure  16.  Particle  source  from  the  floor  as  a  function  of  time.  The  particle 
source  was  calculated  using  the  resuspension  rate  given  by  the  empirical 
correlation,  EC  I.  The  particle  size  was  0.5  (.1111,  and  the  initial  deposition 
amount  was  352  ng  (0.06  ng/cm2).  The  floor  material  was  linoleum. 

For  the  same  condition,  Figure  17  shows  the  temporal  and  spatial  evolu¬ 
tions  of  the  particle  concentration  in  the  plane  of  the  floor  and  in  the  plane 
along  the  diagonal  of  the  chamber  for  4-30  seconds;  Figure  18  shows  the 
same  temporal  and  spatial  evolutions  for  60  -  240  seconds.  Figure  19 
shows  temporal  variations  of  particle  concentration  at  the  exit  of  bottom 
hole  (for  the  APS).  The  particle  concentration  was  predicted  to  reach  a 
maximum  in  a  very  short  time  scale  (~<  60s)  and  then  to  slowly  decrease 
as  the  source  term  became  smaller  (as  the  resuspension  rate  decreased 
over  time)  and  particle  loss  through  the  flow  outlets  continued. 

4.4  Comparisons  of  Resuspension  Data  between  Simulations  and 
Measurement 

Section  4.3.2  (p  48)  presented  the  CFD  simulation  results  of  the  particle 
concentration  for  0.5  pm-sized  particles  over  time.  The  simulations  dem¬ 
onstrated  that  the  empirical  correlation  of  the  particle  resuspension  rate 
could  be  coupled  with  CFD  based  code  as  a  source  term  and  the  temporal 
and  spatial  evolution  of  the  particle  concentration  could  be  obtained  using 
CFD.  In  reality,  all  particles,  whether  biological  warfare  agents  (BWA)  or 
common  dust  particles,  have  size  distributions.  In  the  experiments  on 
which  the  CFD  simulations  were  based,  the  APS  measured  the  particle 
concentration  in  51  different  size  bins  in  the  range  of  0.5  -  20  pm  and  so¬ 
dium  fluorescein  powder  showed  a  size  distribution  between  0.5  and 

15  pm. 
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Figure  17.  Isocontours  of  particle  concentration  in  a  plane  at  the  floor  and  in 
a  plane  along  the  diagonal  of  the  chamber  showing  spatial  and  temporal 
evolution  of  particle  concentration  between  4  and  30  seconds:  at  20  L/min 
(left)  and  132  L/min  (right).  The  particle  size  was  0.5  jim,  and  the  initial 
deposition  amount  was  352  ng  (0.06  ng/cm2).  The  floor  material  was 

linoleum. 
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Figure  18.  Isocontours  of  particle  concentration  in  a  plane  at  the  floor  and  in 
a  plane  along  the  diagonal  of  the  chamber  showing  spatial  and  temporal 
evolution  of  particle  concentration  between  60  and  240  seconds:  at  20  L/min 
(left)  and  132  L/min  (right);  the  particle  size  was  0.5  jim,  and  the  initial 
deposition  amount  was  352  ng  (0.06  ng/cm2);  the  floor  material  was 

linoleum. 
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Time  (s) 

Figure  19.  Time  evolution  of  the  particle  concentration  at  the  outlet  located 
5.25  in.  above  the  floor  at  the  center  between  (A)  o  -  600  s  and  (B)  o  -  3600 
s;  the  particle  size  was  0.5  pm,  and  the  initial  deposition  amount  was  352  ng 
(0.06  ng/cm2);  the  floor  material  was  linoleum. 

Therefore,  to  compare  with  the  experimental  results,  the  CFD  simulations 
in  Section  4.3.2  (p  48)  needed  to  be  carried  out  for  ~  50  different  particle 
sizes.  Although  it  is  certainly  possible,  it  would  take  massive  computa¬ 
tional  power  and  long  computing  time  beyond  the  scope  of  the  Phase  I 
work.  However,  a  simpler  calculation  can  be  done  to  compare  with  the  ex¬ 
perimental  results,  if  it  does  not  account  for  redeposition  due  to  gravita¬ 
tional  settling. 

From  Equation  77,  the  mass  flux  of  particles  at  a  given  size  is: 

0.5787^  12311  0.5989 

Mass  Flux  (dp)  =  8.521  x  10“3  03028P  03f9a  02961  x  InitalMass  Depostion  Density  (dp)  Ecl  79 

*  Pp  Z0  A132 


/r\ 
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Then,  the  mass  resuspended  during  l  hr  for  a  given  particle  size  is: 

j  1.2311  0.5989 

Mass  Resuspended  =  8.521  xlCT3 — - 0296fx  ^M*°5787  x  (Mesh  Area  for  w, ) 

Pv'  Z0  ^132  Mesh  Area 

(•3600  ,  nns 

x  Inital  Mass  Deposition  Density  xj^  t  dt  Eq  80 

Because  the  exponent  of  t  is  >i,  l  is  used  instead  of  o  as  the  lower  limit  of 
the  integration  to  avoid  division  by  o.  It  is  the  same  reason  that  the  resus¬ 
pension  at  time  o  was  assumed  to  equal  zero  and  the  resuspension  started 
at  l  second  in  the  CFD  simulations. 

The  experiments  were  conducted  for  the  resuspension  of  particles  in  the 
size  range  1.7  -  3.8  pm  and  5.0  -  10.3  pm,  separately,  from  the  floor  cov¬ 
ered  with  either  linoleum  or  carpet.  However,  the  results  reported  the  data 
from  the  small  (1.7  -  3.8  pm)  and  the  large  particles  (5.0  -  10.3  pm)  com¬ 
bined.  Figure  20  shows  the  initial  mass  deposition  density  for  linoleum 
and  carpet.  The  small  peaks  between  2  and  3  pm  in  both  cases  were  the 
result  of  combining  the  initial  mass  deposition  densities  of  the  small  and 
the  large  particles. 

For  comparison,  Figures  21-23  show  the  prediction  from  EC  I  against  the 
experimental  data  with  various  fan  settings  for  two  different  floor  materi¬ 
als.  The  fraction  of  mass  resuspended  was  calculated  for  individual  parti¬ 
cle  sizes  (the  ratio  of  the  total  mass  resuspended  during  1  hr  to  the  mass  of 
initial  deposition  on  the  floor  for  a  given  particle  size).  For  the  experimen¬ 
tal  results,  the  mass  resuspended  during  1  hr  was  a  measured  value.  For 
the  predictions  using  EC  1,  the  mass  resuspended  was  calculated  using 
Equation  80.  In  some  predictions,  the  lower  and  upper  limits  of  the  inte¬ 
gration  used  were  different  from  those  specified  in  Equation  80  for  the 
reasons  described  below. 

As  mentioned  earlier,  each  resuspension  experiment  lasted  5  hrs,  1  hr  with 
the  fans  off  and  1  hr  for  each  of  four  fan  speeds  (45,  75, 110,  and  140).  Par¬ 
ticle  resuspension  was  not  observed  in  any  experiment  when  the  fans  were 
off  or  at  speed  45.  Figure  21  shows  the  experimental  results  and  the  pre¬ 
dictions  for  the  fan  speed  of  75.  The  measured  fan  flow  rate  was  20  L/min 
for  the  experiments.  The  red  circles  were  the  prediction  using  EC  I  with 
the  friction  velocity  obtained  from  the  CFD  simulation  at  20  L/min  fan 
flow  rate. 
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Figure  20.  The  initial  mass  deposition  density  as  a  function  of  particle  size 

on  (A)  linoleum  and  (B)  carpet. 
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Figure  21.  The  total  fraction  of  mass  resuspended  during  1  hr  of 
resuspension  experiment  for  (A)  linoleum  and  (B)  carpet.  The  fraction  of 
mass  resuspended  was  calculated  for  a  given  particle  size  as  the  ratio  of  the 
mass  of  particle  resuspension  to  the  mass  of  initial  deposition  on  the  floor. 
The  fan  setting  was  75,  and  the  measured  fan  flow  rate  was  20  L/min.  The 
red  and  the  blue  circles  indicate  predictions  using  EC  I  at  20  L/min  and 
132  L/min,  respectively,  of  the  fan  flow  rate.  The  purple  circles  represent  the 
prediction  using  the  empirical  correlation  given  by  Loosmore  (2003)  using 

20  L/min  as  the  fan  flow  rate. 
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Figure  22.  The  total  fraction  of  mass  resuspended  during  1  hr  of 
resuspension  experiment  for  (A)  linoleum  and  (B)  carpet.  The  fraction  of 
mass  resuspended  was  calculated  for  a  given  particle  size  as  the  ratio  of  the 
mass  of  particle  resuspension  to  the  mass  of  initial  deposition  on  the  floor. 
The  fan  setting  was  110,  and  the  measured  fan  flow  rate  was  40  L/min.  The 
red  and  the  blue  circles  indicate  predictions  using  EC  I  at  40  L/min  fan  flow 
rate,  and  tstart  and  tend  specify  the  lower  limit  and  the  upper  limit, 
respectively,  of  the  integration  in  Eq  80. 
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Figure  23.  The  total  fraction  of  mass  resuspended  during  1  hr  of 
resuspension  experiment  for  (A)  linoleum  and  (B)  carpet.  The  fraction  of 
mass  resuspended  was  calculated  for  a  given  particle  size  as  the  ratio  of  the 
mass  of  particle  resuspension  to  the  mass  of  initial  deposition  on  the  floor. 
The  fan  setting  was  140,  and  the  measured  fan  flow  rate  was  70  L/min.  The 
red,  the  blue,  and  the  purple  circles  indicate  predictions  using  EC  I  at 
60  L/min  of  the  fan  flow  rate,  and  tstart  and  ten a  specify  the  lower  limit  and  the 
upper  limit,  respectively,  of  the  integration  in  Eq  80. 


Since  the  particle  resuspension  did  not  occur  before  setting  the  fan  speed 
at  75,  researchers  assumed  that  the  time  =  o  at  the  beginning  of  the  ex¬ 
periments  and  used  1  and  3600  s  as  the  lower  and  the  upper  limit  of  the 
integration  in  Equation  80.  For  linoleum  floor,  EC  I  (using  20  L/min)  pre¬ 
dicted  reasonably  well  for  particles  between  2  and  4  pm.  However,  the  ex¬ 
perimental  results  for  the  particles  larger  than  4  pm  were  much  lower  than 
the  prediction.  Note  that,  for  the  experiments  with  linoleum  floor,  there 
was  the  discontinuity  in  that  data  between  the  small  and  the  large  particle 
experiments  for  the  fan  speed  75.  This  discontinuity  is  the  reason  for  the 
abrupt  decrease  in  the  fraction  of  mass  resuspended  of  the  experiment  for 
the  particles  larger  than  4  pm  in  Figure  21(A);  also  the  experimental  re¬ 
sults  for  the  large  particles  were  thought  to  be  inaccurate.  For  the  carpet, 
prediction  by  EC  I  (using  20  L/min)  was  slightly  higher  than  the  experi¬ 
ment  for  the  particle  larger  than  2  pm  (Figure  2i[B]).  For  the  particles 
smaller  than  1  pm,  the  experimental  results  were  higher  than  the  predic¬ 
tions  on  both  linoleum  and  carpet. 

During  the  flow  field  simulation,  the  disagreement  between  the  measured 
average  velocity  and  the  predicted  velocity  by  CFD  simulations  in  the 
plane  10  cm  above  the  floor  for  a  specified  fan  flow  rate  were  noted.  At 
20  L/min  fan  flow  rate,  the  average  velocity  from  the  CFD  simulation  was 
lower  than  the  measured  values.  Therefore,  the  friction  velocity  obtained 
from  the  CFD  simulation  using  20  L/min  would  be  lower  than  the  friction 
velocity  inferred  from  the  measured  average  velocity.  In  Section  4.3.2 
(p  48),  the  fan  flow  rate  was  increased  up  to  132  L/min  in  an  effort  to 
match  the  measured  average  velocity  at  20  L/min.  However,  at  20  L/min, 
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the  average  velocity  predicted  from  the  CFD  simulation  using  132  L/min 
was  still  lower  than  the  measured  average  velocity  from  the  experiment. 

Nevertheless,  the  friction  velocity  from  the  simulation  with  132  L/min  was 
used  to  predict  particle  resuspension.  The  purpose  was  to  examine  the  ef¬ 
fect  if  the  fan  flow  rate  measurement  was  erroneous  and  the  average  veloc¬ 
ity  measurement  was  accurate  in  the  experiments.  Because  the  friction  ve¬ 
locities  at  132  L/min  is  higher  than  that  at  20  L/min,  as  one  can  expect, 
the  fraction  resuspended  predicted  using  132  L/min  fan  flow  rate  was 
higher  than  that  predicted  using  20  L/min  (Figure  21).  Since  the  predicted 
values  from  132  L/min  were  much  higher  than  those  from  the  experi¬ 
ments,  it  was  decided  that  all  the  predictions  would  be  calculated  using  the 
friction  velocities  from  CFD  simulations  with  the  fan  flow  rate  that  was  the 
same  as  the  measured  fan  flow  rate,  not  with  the  fan  flow  rate  that  agreed 
with  the  measured  average  velocity. 

The  fraction  of  mass  resuspended  was  calculated  using  Loosmore’s  em¬ 
pirical  correlation.  Figure  21  shows  the  results.  The  predictions  from  EC  I 
and  Loosmore’s  empirical  correlation  were  similar  for  the  0.5  pm-sized 
particles.  However,  for  the  15  pm-sized  particles,  Loosmore’s  empirical 
correlation  predicted  the  fraction  resuspended  two  orders  of  magnitude 
lower  than  EC  I. 

Figure  22  shows  the  experimental  results  and  the  predictions  from  EC  I  at 
the  fan  speed  of  no.  The  measured  fan  flow  rate  was  40  L/min,  and  the 
predictions  used  the  friction  velocities  from  CFD  simulations  with 
40  L/min.  Because  the  experiments  at  the  fan  speed  of  110  were  conducted 
right  after  the  experiments  at  75  fan  speed  for  an  hour,  it  was  not  clear 
how  to  set  the  time  span  for  these  experiments.  If  one  assumes  the  previ¬ 
ous  experiment  at  the  fan  speed  of  75  did  not  have  any  influence  on  the 
experiment  at  the  fan  speed  of  no,  the  time  span  for  the  experiment 
should  be  o  -  1  hr.  If  the  previous  experiment  had  affected  the  experiment 
at  the  fan  speed  of  110,  it  would  have  resulted  in  lowering  the  fraction  re¬ 
suspended  (at  the  fan  speed  of  no)  compared  to  an  experiment  at  the  fan 
speed  of  no  that  is  conducted  independently.  In  that  case,  it  may  be  more 
realistic  to  use  the  time  span  of  1  -  2  hr  to  predict  the  experiment  at  110. 
For  the  linoleum  floor,  using  the  time  span  of  1  -  2  hr  agreed  better  with 
the  experimental  results  than  using  the  time  span  of  o  -  1  hr  for  large  par¬ 
ticles  (>5  pm).  In  contrast,  the  time  span  of  o  -  1  hr  provided  comparable 
results  to  the  experiments  on  the  carpet  for  the  particles  larger  than  2  pm. 
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For  the  particles  smaller  than  l  pm,  EC  I  underestimated  the  resuspended 
fraction  regardless  of  floor  materials  and  time  span  used. 

Figure  23  shows  the  experimental  results  and  the  predictions  from  EC  I  at 
the  fan  speed  of  140.  The  measured  fan  flow  rate  was  70  L/min,  and  the 
predictions  used  the  friction  velocities  from  CFD  simulations  with 
60  L/min.  The  experiments  with  the  fan  speed  of  140  followed  the  experi¬ 
ments  with  the  fan  speed  of  110  so  assigning  the  time  span  for  the  experi¬ 
ments  was  also  a  concern.  Three  different  time  spans  (o  -  1  hr,  1  -  2  hr, 
and  2-3  hr)  were  used  for  the  empirical  correlations.  For  the  linoleum 
floor,  again,  the  result  of  1  -  2  hr  or  2  -  3  hr  agreed  better  with  the  ex¬ 
perimental  results  than  the  result  of  o  -  1  hr  for  large  particles.  On  the 
carpet,  the  time  span  of  o  -  1  hr  gave  a  closer  result  to  the  experiment  than 
1  -  2  or  2  -  3  hr. 

In  conclusion,  EC  I  predicted  the  experiments  reasonably  well  for  particles 
larger  than  3-5  pm  and  underestimated  the  resuspension  fraction  for 
particles  smaller  than  1  pm  in  all  conditions.  Because  the  experiments 
were  carried  out  sequentially  with  increasing  fan  speed,  designating  the 
time  span  of  integration  in  Equation  80  for  the  subsequent  experiments 
was  an  issue.  For  the  linoleum  floor,  time  counting  from  the  first  experi¬ 
ment  gave  better  agreement  with  the  experiments.  This  indicated  that  the 
previous  experiment  may  have  influenced  the  results  of  the  following  ex¬ 
periments.  For  the  carpet,  time  counting  from  the  first  experiment  did  not 
agree  well  with  the  experiments.  Instead,  the  time  span  of  o  -  1  hr  worked 
acceptably  for  the  subsequent  experiments.  For  future  experiments,  differ¬ 
ent  fan  speed  experiments  should  be  carried  out  independently  to  remove 
any  ambiguity  in  assign  the  time  for  resuspension. 

Note  that  the  friction  velocities  used  for  predicting  the  resuspension  frac¬ 
tion  were  from  the  CFD  simulations  assuming  a  smooth  surface  without 
accounting  for  the  surface  roughness.  If  the  surface  roughness  were  in¬ 
cluded  in  the  CFD  simulations  for  the  flow  field  calculation,  the  friction 
velocities  would  have  been  different.  Consequently,  it  would  have  resulted 
in  different  predictions  of  the  resuspended  fraction.  The  prediction  did  not 
include  the  particle  redeposition  due  to  gravitational  settling,  which  may 
have  overestimated  the  resuspension  fraction.  It  is  anticipated  that  these 
issues  will  be  addressed  during  Phase  II  of  this  project. 
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5  Summary,  Conclusions,  and 
Recommendations 

5.1  Summary 

In  Phase  I,  Triton  Systems  developed  empirical  correlations  that  can  pre¬ 
dict  the  particle  resuspension  rate  as  a  function  of  time,  the  friction  veloc¬ 
ity,  the  particle  size,  the  particle  density,  the  surface  roughness,  and  the 
Hamaker  constant,  which  represents  van  der  Waals  interaction  between 
the  particle  and  the  surface.  These  parameters  were  selected  based  on  the 
relevant  mechanisms  for  resuspension  phenomena,  and  the  empirical  cor¬ 
relations  were  expressed  in  dimensionless  numbers.  A  physics-based 
model  using  the  moment  balance  between  adhesion  and  removal  forces 
was  also  developed.  The  empirical  correlations  were  tested  against  ex¬ 
perimental  data  available  in  literature.  Finally,  the  CFD  simulations  of  the 
flow  field  were  performed  in  a  chamber  where  particle  resuspension  ex¬ 
periments  were  conducted.  The  CFD  simulations  of  the  particle  resuspen¬ 
sion  and  transport  in  the  chamber  were  also  carried  out  for  one  particle 
size  with  different  floor  materials.  The  resuspended  fractions  predicted  by 
the  empirical  correlation  EC  I  using  the  friction  velocities  from  the  flow 
field  simulations  was  compared  to  the  experimental  results,  and  the  pre¬ 
diction  agreed  reasonably  well  with  the  experiments. 

5.2  Conclusions 

Among  the  seven  empirical  correlations  derived  in  Phase  I,  three  correla¬ 
tions  (EC  I,  EC  IV,  and  EC  VII)  were  physically  realistic.  EC  I  and  EC  IV 
expressed  the  resuspension  rate  as  a  function  of  time,  the  friction  velocity, 
the  diameter  and  the  density  of  particles,  the  surface  roughness,  the 
Hamaker  constant,  and  the  density  of  air.  The  difference  between  EC  I  and 
EC  IV  was  the  method  used  to  assign  the  reference  time  for  the  resuspen¬ 
sion  rate  that  was  calculated  from  the  fraction  removal  over  time.  EC  I 
used  the  mid-point  time  as  the  representative  time,  and  EC  IV  used  the 
time  t'  that  divided  the  area  under  the  l/t  curve  for  a  given  time  interval 
into  two  equal  areas.  The  resuspension  rate  in  EC  VII  was  a  function  of 
time,  friction  velocity,  diameter  and  density  of  particles,  surface  rough¬ 
ness,  and  the  density  and  viscosity  of  air.  EC  VII  employed  the  same  ap¬ 
proach  to  obtain  the  representative  time  for  the  resuspension  rate  as  EC 
IV.  These  three  empirical  correlations  predicted  experimental  data  from 
the  literature  reasonably  well.  For  some  sets  of  data,  the  ratio  of  the  pre- 
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dieted  to  the  observed  resuspension  rates  was  centered  at  l  and  scattered 
within  the  range  of  one  order  of  magnitude  in  either  direction.  However, 
the  correlations  underestimated  the  resuspension  rate  for  other  sets.  The 
predictions  from  the  three  empirical  correlations  were  comparable,  and 
their  performances  were  acceptable  overall. 

The  models  that  used  t'  as  the  representative  time  were  expected  to  yield 
better  correlations  than  those  that  used  the  mid-point  time.  In  fact,  using 
t'  did  not  show  any  advantageous  results  over  using  the  mid-point.  Al¬ 
though  EC  I,  EC  IV,  and  EC  VII  are  equally  valid  to  predict  the  resuspen¬ 
sion  rate,  researchers  slightly  preferred  EC  I  due  to  the  fact  that  it  uses  the 
simple  mid-point  for  the  representative  time  and  includes  the  Hamaker 
constant,  which  accounts  for  the  material  interactions  between  a  surface 
and  particles.  One  disadvantage  of  EC  I  is  that  the  Hamaker  constants  for 
some  materials  are  not  readily  available. 

5.3  Recommendations 

The  quality  of  an  empirical  correlation  depends  on  the  data  used  to  draw 
the  correlations.  Therefore,  to  model  and  predict  the  indoor  resuspension 
accurately  it  would  be  essential  to  derive  models  based  on  the  experimen¬ 
tal  data  with  the  conditions  that  are  characteristic  to  indoor  resuspension. 
The  empirical  correlations  developed  in  Phase  I  can  be  refined  with  better 
experimental  data. 

The  physics-based  model  suggests  the  group  of  dimensionless  numbers 
that  were  slightly  different  than  those  obtained  for  the  empirical  correla¬ 
tions.  It  is  worthwhile  to  derive  another  empirical  correlation  based  on  the 
dimensionless  numbers  from  the  physics-based  model  and  compare  the 
performance  with  the  empirical  correlations  obtained  in  Phase  I. 
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Acronyms,  Abbreviations,  and  Terms 


Term 

AFM 

APS 

BWA 

CB 

CERL 

CFD 

CFDRC 

DC 

EC 

EC 

ERDC 

ESI 

FC 

FD 

FY 

HVAC 

IAQ 

JKR 

LBNL 

MIPR 

NCRP 

NRPB 

PVC 

RH 

RH 

SBIR 

SS 

TNT 

TR 

UK 

URL 

USAEC 

WWW 


Spellout 

atomic  force  microscopy 

Aerodynamic  Particle  Sizer 

biological  warfare  agents 

chemical  and  biological 

Construction  Engineering  Research  Laboratory 

computational  fluid  dynamics 

CFD  Research  Corporation 

direct  current 

Engineer  Circular 

Empirical  correlation 

Engineer  Research  and  Development  Center 
Environmental  Sustainability  Index 
Facility  Composer 
Forced  Draft 
fiscal  year 

heating,  ventilating,  and  air-conditioning  (FIVAC) 
indoor  air  quality 

Johnson-Kandall-Roberts  (model) 

Lawrence  Berkeley  National  Laboratory 

Military  Interdepartmental  Purchase  Request 

National  Council  on  Radiation  Protection  &  Measurements 

National  Radiological  Protection  Board 

polyvinyl  chloride 

relative  humidity 

Relative  humidity 

Small  Business  Innovative  Research 

Stainless  Steel 

trinitrotoluene 

Technical  Report 

United  Kingdom 

Universal  Resource  Locator 

U.S.  Army  Environmental  Command 

World  Wide  Web 


A 132  Hamaker  constant  for  Particle  1  and  Surface  2  in  Medium  3 
a  Rate  constant 

a 0  Coefficient  in  empirical  correlation 
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Term  Spellout 

a  1  Coefficient  in  empirical  correlation 

a2  Coefficient  in  empirical  correlation 

bO  Coefficient  in  empirical  correlation  in  EC  I  and  EC  IV 

dl  Coefficient  in  empirical  correlation  in  EC  I  and  EC  IV 

b2  Coefficient  in  empirical  correlation  in  EC  I  and  EC  IV 

b3  Coefficient  in  empirical  correlation  in  EC  I  and  EC  IV 

M  Coefficient  in  empirical  correlation  in  EC  I  and  EC  IV 

cO  coefficient  in  empirical  correlation  in  EC  II  and  EC  V 

cl  Coefficient  in  empirical  correlation  in  EC  II  and  EC  V 

c2  Coefficient  in  empirical  correlation  in  EC  II  and  EC  V 

c3  Coefficient  in  empirical  correlation  in  EC  II  and  EC  V 

CfX  Local  shear  stress  coefficient 

c/0  Coefficient  in  empirical  correlation  in  EC  III,  EC  VI,  and  EC  VII 

c/1  Coefficient  in  empirical  correlation  in  EC  III,  EC  VI,  and  EC  VII 

c/2  Coefficient  in  empirical  correlation  in  EC  III,  EC  VI,  and  EC  VII 

c/3  Coefficient  in  empirical  correlation  in  EC  III,  EC  VI,  and  EC  VII 

c/4  Coefficient  in  empirical  correlation  in  EC  III,  EC  VI,  and  EC  VII 

c/5  Coefficient  in  empirical  correlation  in  EC  III,  EC  VI,  and  EC  VII 

c/h  Hydraulic  diameter 

c/p  Particle  diameter 

Ei  Young’s  modulus  for  particle 

E2  Young’s  modulus  for  surface 

F  Ratio  of  adhesion  force  to  removal  force 

Eadh  Adhesion  force 

Ec  Capillary  force 

Ed  Drag  force 

Eedl  Electrostatic  force  due  to  electrostatic  double  layer 
El  Lift  force 

Fq  Electrostatic  force  due  to  a  charged  particle 

Evdw  van  der  Waals  force 

g  Gravity 

K  Resuspension  factor 

k  Turbulent  kinetic  energy 

M  Ratio  of  adhesion  moment  to  removal  moment 
N  Particle  number  density  on  surface 

No  Number  density  of  particles  on  surface  at  time  0 
Nt  Total  particle  number  density  on  surface 

pt  Turbulent  pressure 

q  Charge  on  particle 

ra  Contact  radius 


ra 0  Contact  radius  at  zero  applied  load 

Adhesion  Adhesion  radius 


ERDC/CERL  TR-08-4 


65 


Term 


h 

s 

t 

Oo 

u* 

Vd 

Zo 

£ 

so 

(f)c 

V 

yiv 

V 

K 

A 
M a 

V 
Vl 
V2 
Pa 

Pp 

<y 

Tvj 


Spellout 
Particle  radius 

Separation  distance  between  the  particle  and  the  surface 
Time 

Free  stream  velocity 

Friction  velocity 

Turbulent  deposition  velocity 

Aerodynamic  surface  roughness  length 

Turbulent  intensity 

Dielectric  constant 

Contact  potential 

Shear  rate 

Liquid-gas  surface  tension 
Diabatic  influence  function 
von  Karman  constant 
Resuspension  rate 
Viscosity  of  air 
Kinematic  viscosity  of  air 
Poisson  ratio  of  particle 
Poisson  ratio  of  surface 
Density  of  air 
Particle  density 
Surface  energy 
Wall  shear  stress 
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